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TRACKING OF THE ATS-3 SYNCHRONOUS SATELLITE 

BY THE 

VERY LONG BASELINE INTERFEROMETER (VLBI) TECHNIQUE 
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Geodynamics Branch 

Trajector Analysis and Geodynamics Division 
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SUMMARY 

During 1971, a series of very long baseline interferometer 
(VLBI) observations were made of the C-band (6 cm) radio 
signals from the ATS-3 communications satellite which is 
in a synchronous, near- equatorial orbit. The first series 
of observations were conducted during May-June 1971 from 
Rosman, North Carolina (NASA/ ATS Station 85' dish) and 
Mojave, California (NASA/ ATS Station, 40' dish). The second 
series of observations were conducted during August- 
September, 1971 from Rosman, North Carolina (NASA/ ATS 
Station, 85' dish), Owens Valley, California (Cal Tech, 130' 
dish) and Agassiz, Massachusetts (SAO Agassiz Radio 
Observatory, 84' dish). The ATS-3 Spacecraft position 
has been determined with a precision of 70-100 meters 
and its velocity with a precision of less than a mm/ sec. 

The ATS-3 orbital elements have been computed using the 
GEODYN program and the derived values are consistent 
with those derived from conventional tracking data. 
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1 . 


INTRODUCTION 


During 1971, a series of very long baseline interferometer (VLBI) 
observations were made of the C-band (6 cm) radio signals from the ATS-3 
communications satellite which is in a synchronous, near-equatorial orbit. The 
first series of observations were conducted during May-June 1971 from Rosman, 
North Carolina (NASA/ATS Station, 85' dish) and Mojave, California (NASA/ATS 
Station, 40' dish). The second series of observations were conducted during 
August - September, 1971 from Rosman, North Carolina (NASA/ATS Station, 

85' dish), Owens Valley, California (Cal Tech, 130' dish) and Agassiz, 
Massachusetts (SAO Agassiz Radio Observatory, 84' dish). The ATS-3 Space- 
craft position has been determined with a precision of 70-100 meters and its 
velocity with a precision of less than a mm/sec. The ATS-3 orbital elements 
have been computed using the GEODYN program and the derived values are 
consistent with those derived from conventional tracking data. 
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2. EXPERIMENT DESCRIPTION 

2.1 EXPERIMENT #1 

During May - June, 1971, the first series of interferometric observations 
were made of the C-band (6 cm) radio signals from the ATS-3 communications 
satellite which is in a synchronous, near equatorial orbit. Figure 1 shows the 
schematic set-up of the experiment in which Rosman and Mojave ATS stations 
were used as an interferometer pair. Some of the experimental parameters 
are given in Table 1. 


TABLE 1 


ATSR Antenna Parameters - Experiment # 1 


System Parameter 

Rosman 

Mojave 

Diameter 

85' 

40' 

Type of Mount 

X-Y 

X-Y 

System Noise Temperature 
(Cooled Parametric Amplifiers) 

65° K 

65° K 

Receiver Gain 

58 db 

51 db 

Receiver Beam Width 

0.2° 

0.47° 

Signal Frequency 

4178.59072 MHz 

4178.59072 MHz 

rs ... ( Latitude 

eo e ic ) Longitude 

Coordinates ) Altitude 

35° IT 56" .07 N 
277° 7' 27" .19 E 
828 meters 

35° 19' 53" .78 N 
243° 6' 43" .86 E 
887 meters 


Scheduling at the tracking stations permitted normal observational runs 
of up to four hours per day. On several occasions, extended observations were 
made for over 12 hours. Most of the observations were with the ATS-3 satellite 
even though a few observations were made with the ATS-1 satellite. 

There are two C-band transponders on the ATS-1 and ATS-3 satellites. 
Transponder #1 receives at 6212.094 MHz and transmits at 4178.59072 MHz. 
Transponder #2 receives at 6301. 050 MHz and transmits at 4119. 5990 MHz. It is 
transponder #1 that was used in most of our observations. The orbital para- 
meters of the ATS-1 and 3 satellites are given in Table 2. 
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Figure I . Schematic Set-up of VLBI Experiment 
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TABLE 2 


ATS-1 & ATS-3 Orbital Elements 
(Furnished by the GSFC ATS Projects Office) 


Orbital Element 

ATS-1 

ATS-3 

Apogee (km) 


35813.15 

35910.12 

Perigee (km) 


35760.045 

35664.40 

Longitude (degrees) 


149° . 152 

78°. 850 

Latitude (degrees) 


0° . 02 1 

1° . 103 

Period (minutes) 


1436.101 

1436.018 

Drift Rate (deg/day) 


0° . 008 E/day 

0°.013E/day 

Inclination (deg) 


3°. 349 

1 ° . 706 

Eccentricity 


0.00625 

0.00291 

Semi-major Axis (Km) 


42164.96 

42165.43 

Spacecraft Attitude | 

Right Ascension 
Decl ination 

161°. 5 
85°. 6 

165°. 5 
88°. 3 


The satellite transponder radiated signals with an EIRP (effective iso- 
tropic radiated power) of 446 watts centered at 4178.59072 MHz and extending 
over a 20 MHz bandwidth. The signal is so strong - nearly six orders of mag- 
nitude larger than those from natural celestial sources, one could have used 
very small (~3’) microwave horns in the place of the large parabolic antennas 
that are available at the sites. The quality of the satellite fringes is, as expected, 
superior to those from the quasars due to the enhanced signal-to-noise ratio. 

2.2 EXPERIMENT #2 

The second experiment was performed during August - September 1971. 

In this experiment, simultaneous observations were made from three stations: 
Rosman, North Carolina (ATS-85' dish); Owens Valley, California (Cal Tech- 
130’ dish); Agassiz, Massachusetts (SAO-84' dish). The quasar observations 
were made at 4995 MHz and the ATS-3 satellite observations were made at 
4178 MHz. Results of the second experiment are reported separately. 
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2.3 THE INTERFEROMETER SYSTEM: 

A VLBI interferometer system consists of a pair of terminals (Rosman 
and Mojave in this case) with the following elements (Figure 2). 

• Antenna and feed 

• Radiometer/receiver 

• IF interface and calibrating equipment 

• VLBI Backend (data conditioning and recording equipment) 

• Timing and Frequency control system 

• Correlator/ Computer 

2.3.1 Antenna/Receiver System: 

Antenna and receiver specifications at Mojave and Rosman have been 
discussed earlier. 

2.3.2 IF Interface and Calibrating Equipment 

The receiver IF frequency at the Rosman and Mojave stations is 70 MIIz. 
However, the video converter of the VLBI backend is designed to accept either 
a 30 MHz signal or a 150 MHz signal as an input. In order to derive the required 
30 MHz signal, suitable "IF Interface" equipment is used. This interface equip- 
ment is merely a synthesizer/mixer combination which uses a 40 MHz signal to 
mix with the station IF of 70 MHz and thereby generates a new 30 MHz IF (Figure 3). 

2.3.3 VLBI Backend and Recording System 

The principal components of the instrumentation system are the video 
converter, the sampler formatter and the video recorder. Figure 4 is a diagram 
of the instrumentation system. One rack contains the time-keeping equipment 
and the frequency source and the other contains the video converter and the 
data-formatting circuitry. The latter is modularly assembled from DTL-logic 
integrated -circuit cards. 

The backend system accepts either a 30- or 150-MHz IF signal from the 
receiver. After two or three stages of phase-coherent conversion, the filtered 
..■ video signal is amplified, limited, and clipped to provide only zero crossing in- 
formation within the selected video passband. 
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Figure 2. Block Diagram of an Interferometer System 
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Figure 3. IF Interface Equipment 


The nonreturn -to-zero (NRZ) sampler samples 40-nsec windows of the 
clipped video at a rate controlled by the frequency standard which also provides 
the basic frequency for phase coherent conversion of the IF to video frequency. 
After being sampled, the 1-bit data (zero if sampled below the baseline, one if 
sampled above) are serially applied to a multiplex circuit that arranges the 
format of the data frame. At the start of each video frame, the multiplex unit 
inserts a recorder head-switching period (all ones), a 60 bit synchronization 
word (all zeroes), and a BCD time-code word. These nondata bits are generat- 
ed at the same sampling rate as are the data bits, and their position within the 
frame is fixed by the clock counter and the event decoder. The time-code word 
is transferred once each frame from the digital clock into a register and is then 
transferred to the multiplexer. A recording "start time" is preset by thumb- 
wheel switches. Recording is automatically initiated upon receipt of the first 
frame time occurring after the selected "second" of time. 

The data frame is encoded in the modulator into a self-clocking 
recording format. Either bi-phase return to zero (RZ) or Manchester Coding 
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Figure 4. The Instrumentation System 
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may be selected. The videotape recorder used is a modified Ampex 660-C 
helical scan TV-type recorder with provision for external head-drum 
synchronization. 

2. 3. 3.1 Video Converter, Clipper/Sampler. The basic circuitry of the video 
converter is shown in Figure 5. 


WIDEBAND AMPLIFIER 



5 MHz INPUT 
FROM FREQUENCY 
STANDARD 


1 MHz INPUT 


Figure 5. Video-Converter 

The intermediate frequency input to the video converter is single side- 
band filtered around a center frequency of either 30 or 150 MHz. The video 
converter converts the IF from the low-noise receiver into a video signal 
suitable for 1-bit digital sampling and subsequent recording. The IF signal is 
amplified and filtered; the phase is coherently converted in frequency to estab- 
lish a single side-band output signal in which the carrier frequency is zero (dc). 
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The output video signal band extends from a lower frequency cut off of 50 kHz 
to an upper frequency (3-db) cut off of one of the four selectable steps: 350 kHz, 
480 kHz, 960 kHz and 2.0 MHz. Uniform response is achieved throughout the 
band with sharp rejection at the edges using appropriate filtering process. 

The video converter contains the necessary frequency multiplier, fre- 
quency synthesizers and phase -locked oscillator circuitry to generate 120, 27 
and 2 MHz phase -coherent local oscillator signals from the 1 and 5 MHz signals 
supplied by the precision frequency standard. The 150 MHz signal is amplified 
by a broad band (100-200 MHz) amplifier and fed to a 150/30 MHz mixer. 

When the video converter is operated with a 30 MHz input signal, the IF 
signal is fed through a broad band amplifier that acts as a gain-control element 
in the video converter. The 30 MHz signal is fed to a 30/3 MHz mixer. The 
local-oscillator signal at 27 MHz is obtained from a phase-locked oscillator. 

The reference or phase-locking signal for this 27 MHz oscillator is obtained by 
two multipliers (x3 and x9) in sequence driven by a 1 MHz reference signal from 
the frequency standard. An unlock oscillator circuit is connected to the phase- 
locked oscillator to provide automatic frequency searching and visual indication 
of loss of phase -lock. 

The output band from the 30/3 MHz mixer is fed to a series of bandpass 
amplifiers centered on 3.0 MHz. Each bandpass amplifier utilizes thyee tuned 
LC circuits to achieve a uniform response in the 2. 0-4.0 MHz band with high 
rejection below 2.0 MHz and above 4.0 MHz. Tuned LC trap circuits between the 
first and second and the second and third bandpass amplifiers provide additional 
rejection in the frequency region just below 2.0 MHz. 

The output of the band pass amplifier is fed to the 3/1 MHz mixer, where 
a 2 MHz local oscillator signal is mixed with the 2 to 4 MHz signal band. The 
2.0 MHz local oscillator signal is derived from a frequency doubler driven by 
the 1 MHz signal from the frequency standard. A video amplifier with a lower 
frequency cut off of approximately 25 kHz and an upper frequency cut off of 
approximately 1.8 MHz is a part of the mixer-output circuitry. The output of 
the video amplifier/filter is a single sideband signal that has a carrier frequency 
of zero. An increase in frequency in the IF (30 or 150 MHz) will provide a 
corresponding increase in the frequency of the video band. 

The diode limiter removes a major portion of the amplitude fluctuation 
in the video signal. The broad band noise input produces the same signal level 
on each of the bandwidth settings through internal adjustments of the relative 

gain. 
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The signal output from the limiter is connected to the clipper circuitry. 
Here the video-band signal is amplified and clipped in amplitude so that only the 
zero-crossing information is retained. For a noise input signal, the output from 
the clipper is a random square wave in which the voltage level is either zero or 
+ 4 volts. 


2.3. 3.2 Data Recording System. The data format used in recording the 1-bit 
digital data on videotape is illustrated in Figure 6. The individual data bits 
consist of self-clocking bi -polar or phase-modulated pulses, depending upon the 
modulation mode selected. 


12,000 BITS - ONE VIDEO BLOCK OF DATA 


0-23 

1 

1 

24-83 

84-137 

« [ 

i 

138-11,975, ! 


(head- 

(all zero 

( time- 

(data bits ) 

! 


switching 

sync word ) ’ 

code ) 



ALL ONES 


11,976-11,999 


EXTERNAL 

TIME-CODE 

REGISTER 

LOCKS-UP 


RECORD SYSTEM 
TIME-CODE REGISTER 
LOADS 


DATA BITS 
BEGIN 


a) EXTERNAL TIME CODE 
REGISTER UNLOCKS 

b) RECORD SYSTEM 
TIME-CODE REGISTER 
RESETS 


Figure 6. Video-Tape Format 

A frame of data occupies the time required for one pass of the record 
head across the tape, or 16-2/3 msec. At a nominal sampling rate of 720 
kbits/sec, a frame would accomodate 12,000 bits. Of these, 11,838 bits represent 
data; the remaining 162 are used for time-code and housekeeping words. At 
the beginning of each frame, the vertical sync signal corresponds to the begin- 
ning of a 24 bit all -one sequence designed to bracket the head -switching transient 
of the video play-back. Thus transients occurring during this period do not 
affect the data. Following the 24-bit sequence, a 60-bit all zero sequence is 
written. This serves as an unambiguous synchronization and frame-identifier 
during playback and is used to control the reformatting on the digital tape. 
Immediately after the 60-bit sequence, a 54-bit 8-4-2-1 BCD time-code, which 
registers the time of occurrence of the start of the frame (coinciding with the 
vertical sync pulse), is recorded. The time code extends from day of year to 
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units of microseconds. The data bits follow the time code. Prior to the end of 
the frame, 24 more ones are inserted to complete the head-switch word that 
brackets the start of the frame. Because each pulse is precisely related to each 
other pulse and because the time-code tags the start of the sync-word, recorder 
time-base errors in play-back (or recording) are not important. Also, signals 
to lock and unlock the external clock and load the recording system's time-code 
register are generated coinciding with the boundaries of the words in the frame. 

The circuitry required to record the format described above on a video 
recorder is presented in Figure 7. 



WORD SWITCHING 


Figure 7. Data-Formatter Circuitry 

The sampling frequency (720,960 kbits/sec) is generated either by 
direct division of the 5 MHz standard frequency or with the use of a nonintegral 
divider or frequency synthesizer. This sampling signal is used both to toggle 
the NRZ sampler flip-flop (which provides the clipped data to the multiplexer) 
and to drive the sample clock counter. Outputs of this counter are selected and 
combined to produce a number of logical commands. These commands trigger 
the multiplexer gates that assemble the various words that comprise the data 
frame. The decoding process which is all done on a bit-count basis, produces 
consecutively the trigger signals for the vertical sync, first head-switch, 
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sync-word, time-code, data-bits, and the second head-switch. At the end of the 
frame, a reset signal is initiated. 

All bits that reach the modulator are generated in the multiplexer by 
gated flip-flops. When the system switch position is "Stop", the modulator 
receives all ones. When the switch position is "GO", the input to the modulator 
is all ones until the "Start time" trigger pulse generated by the coincidence of 
the clock reading with the preset "Start time" sets a flip-flop. At this time, the 
sync word is written on the tape. It is not possible to generate a "Start time" 
trigger pulse when the system switch is in the "Stop" position. 

The modulator generates both bipolar RZ and phase -modulated (Manches- 
ter Coding) self -clocking record signals and feeds either type signal to the 
record heads, depending on the choice made by the operator. Thp present bit- 
rates are considerably lower than the maximum that can be handled by the 
recorder. Therefore, the bipolar output is generally used since it can be 
visually decoded. At bit rates above 1 Mbit/sec, the phase-modulated output 
will be advantageous. For the bipolar output, the NRZ data are clocked into 
two separate RZ signals: one with pulses representing data ones and the other 
with pulses representing data zeroes. 

The video recorder is an Ampex VR-660C model, dual-head, helical scan 
machine with 2" tape width space designed for color television use. The record- 
er has large bandwidth, high recording density, and excellent servo control. 

These features make it possible to use the 66 OC to record and play back high- 
speed digital data economically on a machine designed for analog data. Servo 
control of the head drum, Capstan drive, and tape tension permits the data frames 
to be precisely synchronized with one head -switching rate of the recorder during 
recording and playback. The recorder, which has a 4.2 MHz bandwidth, can 
record conventional digital data at rates up to 1.5 Mbit/sec. The data rates can 
be as high as 5 Mbit/sec with the use of special encoding techniques (for 
example, sinusoidal wave-phase encoding). Recordings up to 4 hours duration 
can be made on a single reel of tape which is 4,800 feet long. 

In the recording process, timing information derived from the system 
clock is used to automatically start a recoding program and also to incorporate 
in each data frame a 54 bit time word with a precision of 1 ft sec. For the 
first step, thumbwheel switches are used to set the desired start time in terms 
of hours, minutes and seconds. Recording of data begins when the clock reading 
coincides with this setting. Since the clock is interrogated at the frame rate, 
the first frame will appear in less than 16.67 msec after the preset recording 
time. Starting with the first frame, the time word is inserted. The units of time, 
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ranging from hours to microseconds, are transferred from the clock register 
to the time code register and then to the tape. 

2.3.3. 3 Data Playback System. The playback system demodulates the data 
frame and places the bits in a shift register prior to shifting them out to the in- 
puts of a nine-track digital recorder. The play-back instrumentation is shown 
in Figure 8. - 



Figure 8. Playback Instrumentation 


The arranging and shifting of the demodulated bits are predicted on the 
detection of the sync word. Similarly, the initial and terminal frames of data 
transcribed onto the digital tape are determined from a comparison of each 
time-code word with preset digital switch settings. Various formatting and 
housekeeping functions for the digital recording are derived from the clock 
counter and event decoder, which also provide the strobe signals to the recorder. 

The digital data frame is transcribed from the video recorder track on 
to a nine-track digital tape, using a standard (ASC2) computer data format. The 
digital tape format is shown in Figure 9. 
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ONE DIGITAL RECORD= 12 VIDEO FRAMES 


Figure 9. Digital Tape Format 

Since the recorder has a maximum bit-packing density of 800 bits per 
inch and a tape speed of 150 inches/sec, the maximum data bit rate per track 
is 120 kbits/sec. 

Therefore, the maximum data rate for the 9-track recorder is 960 kbit/ 
sec; one track is used for vertical parity. The system is designed to accomo- 
date the lower data-sampling rate of 720 kbit/sec. For data bandwidths greater 
than 960 kbits/sec, the bit packing density would have to be increased. 

A digital-tape data frame consists of 2000 bytes (1 byte = 6 bits) 
recorded during 16-2/3 msec. It begins with a 10-byte sync word followed by 
a 9-byte time-code, 1973 bytes of data and an 8-byte word to span the period 
of the head-switching transients at the end of the frame. Time-base error in 
the video playback can be corrected by varying the length of the last word so 
that the digital-data rate of 120 kbyte/sec is satisfied. Each digital data frame 
corresponds to a video data frame. One digital record is composed of 12 data 
frames and is therefore 200 msec long. Part of every thirteenth video-data 
frame is omitted so that a 5 msec interrecord gap can be written. 


The bit by bit operation of the playback system is self-clocked. The clock 
frequency is derived from the output of the demodulator, with the exception of 
the head-switch transient period, during which time a phase-locked oscillator 
is used to maintain the data-clocking frequency. Besides clocking the data 
through the 60 bit register and driving the counter, the clock output is also 
divided down to the 120 kbit/sec rate to drive the write strobe of the digital 
recorder. 
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Various tape interface functions such as end-of-file (EOF), vertical 
parity check, and check -word characters are also generated by the playback 
circuits. Provision is also made to reset manually the decoding sequence if the 
synchronization word should fail to be detected subsequent to the start of play- 
back. The digital recorder is a 9-track Ampex TM-12 high speed recorder. In 
its record mode, all control functions are assumed by the playback-system 
electronics. 

2. 3. 3. 4 Frequency and Timing System. The sampling and recording of the data 
on the video recorder are controlled by a precision frequency standard which 
also provides the frequencies used in the coherent conversion of the radio 
frequency input signal to video. Either a Hydrogen maser or a rubidium vapor 
cell (HP-5065A) oscillator is used as the frequency standard. At C-band, the 
phase stability is less than one radian in 100 seconds with the rubidium standard. 
This is improved by almost a factor of 50 with the maser. While the maser 
is very essential for precision measurements of fringe frequency, the rubidium 
oscillator is sufficiently stable for the levels of time and phase- re solution 
currently required. 

The frequency standard must also drive clocks and maintain a time- 
synchronization between stations of at least 5 p sec during the experiment. 

Figure 10 shows the VLBI timing system. 

The frequency standard provides a 1 MHz drive to the digital clock, which 
accumulates and visually displays the Universal time with 1- fi sec precision. 
Associated with the clock is a delay generator that accepts the system 1 pps 
signal and digitally delays it up to 999,999 /xsec. The purpose of the delay 
generator is to set in the Loran-C propagation delay associated with various 
Loran stations. It can also be used to set other time signals such as the 1 pps 
from a portable clock in coincidence with the system 1 pps, thereby determining 
their relative offset. 

The principal clock-setting means are comparisons with portable clocks 
carrying accurate UTC (USNO) time and with Loran-C transmissions. The 
digital clock is designed so that it can be directly reset to the 1 pps portable 
clock signal. Alternatively, the portable clock can be used to demand a reading 
of the system clock at a 1-pps rate, thereby displaying the epoch difference 
between the two clocks. When it is desired to set the system clock with Loran-C, 
a selected zero-crossing of the received radio wave form is displayed on an 
oscilloscope. The scope trigger is provided by the Loran Rate Generator and is 
the simulated Loran pulse rate delayed by the propagation time. The delayed 
pulse rate is produced by resetting the rate generator with a pulse generated at 
a time of coincidence (TOC), suitably delayed to account for propagation and 
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system delays. The TOC is determined from Naval Observatory tables for a 
particular Loran chain. On the east coast, it is Cape Fear, North Carolina. 
On the west coast, it is Dana Point, Indiana. After this resetting, the rate- 
generator output is coincidental with the time-of arrival of the -Loran-C pulses. 
After visual alignment is attained with the oscilloscope, the Loran 1 pps is 
used either to reset the clock or to determine the epoch difference. 



Figure 10. VLBI Timing System 


The "demand" inputs of the system clock permit the clock display to be 
frozen at receipt of the demand input pulse. Provision is made to "demand" the 
clock with the Loran - generator output, the delayed clock 1 pps, a VLF 1 pps 
signal and the VLBI time-demand signal from the event decoder. In addition to 
a Loran receiver, a VLF receiver and chart recorder are included in the timing 
system to provide a means of monitoring long-term frequency drift of the 
frequency standard. 

In the Rosman-Mojave experiment, the VLBI clock at Mojave was set by 
a portable Cesium clock and was also regularly monitored by a microwave 
relay signal, both derived from neighboring JPL/Goldstone Echo site. The 
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Rosman VLB clock was set by Loran-C transmission from Cape Fear, North 
Carolina. 

The two VLBI clocks were also synchronized using the satellite pulse 
exchange technique. Pulses at C-band with very sharp rise-time and of 10 
H sec duration were exchanged by the two stations through the dual trans- 
ponders of the ATS-3 geostationary satellite. At each station, a time-interval 
counter was started by the transmitted pulse and stopped by the pulse received 
via satellite from the other station. The probable error of the clock-offset as 
measured by the counter is 10 nanoseconds. The clock offset data was recorded 
regularly. Whenever the offset value exceeded 5 p sec, the Rosman clock was 
adjusted using the delay -generator to keep the offset well within 5 ^ sec. This 
technique worked out very well and added great reliability to the experiment. 



ATS-3 

SATELLITE 



Figure 11. Geometrical Configuration of a Single 
Baseline Experiment 
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3. INTERFEROMETRIC EQUATIONS 


In the simplest terms, an interferometor is an instrument for measuring 
the phase difference between electromagnetic waves received simultaneously at 
any two of its terminals. The very long baseline interferometer (VLBI) uses 
atomic frequency standards to preserve electrical phase coherence between its 
terminals, thus dispensing with the cables or microwave links customarily em- 
ployed for this purpose. Closely regulated, but free-running and independent, 
these standards can be placed at arbitrarily chosen locations to allow the base- 
lines to achieve a length of the order of several million wavelengths of the 
received signal. Given such long baselines, low-noise front-end receivers, wide- 
band VLBI recording systems and stable local oscillators, the VLBI systems are 
able to resolve signal time-delays to an accuracy that is several orders of mag- 
nitude better than that is accomplished by conventional tracking techniques. 

The VLBI technique has been demonstrated to have the capability to 
resolve radio sources with a precision of the order of a millisecond of arc. 

The potential for precision measurement of angular positions of radio sources 
is correspondingly great and is just now being exploited. The adaptation of 
VLBI to satellite tracking therefore appears to be a highly desirable objective. 

Figure 11 shows the geometrical configuration of a single baseline 
VLBI experiment. From the geometry 


d + R-2 = R-2 cos 7 + D cos 9 

(la) 

R2 sin 7 = D cos 9 

(lb) 


The difference in the time of reception of the signal at the two antennas 
is then given by 


or 


£ 

c 



cos 9 - 


( 1 - cos 7) 
sin 7 


sin 9 


(2) 


D cos (9 + 7/2) 
g c cos 7/2 


( 3 ) 


Preceding page blank 
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since 



1 - cos 7 
sin 7 


and c denoting the velocity of light. 

If D is known and r g and R 2 are measured then by Eqs. lb and 3, the 

angle 9 is determined. When the satellite recedes to infinity, 7^0 , the time 
0 D cos 0 

delay expression takes the form applicable to astronomical sources, Tg = — 

An alternative for the satellite time delay is 


r g 


_ R i 


R, 


( 4 ) 


where R[ and R 2 are station to satellite range for the respective stations. 
The fringe frequency is given by 


v t = f 


dt 


( 5 ) 


‘where f is the signal frequency this can be written as a relative doppler 
frequency 

"f = 7 ”^ 2 ) (6) 

Since we are interested in the nature of the correlation function, it is 
necessary to analyze the interferometer response to wideband noise. 


The voltage induced by the source in the antenna #2 is of the form 


V 2 (t) = e i l w(t - T E> + 0T2 ^] 


( 7 ) 


Where u> is the received signal frequency and 0 T2 is the phase shift due to the 
•transmission media linking station #2 and the source. r g is the time-delay 
or the differential time of arrival of the identical wavefront at station #1 relative 
to the arrival time at Station #2. 
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The voltage induced by the source at the antenna #1 is of the form 


Vj(t) = e i [ a,t + 0 Ti(‘)] 


(8) 


The above expressions signify that the wavefront reaching station #1 at a time 
t has reached station #2 at an earlier time, (t -r ). 0 T1 (t) is the phase shift 
due to the transmission media linking station #1 and the source. 

The received signals are then mixed with corresponding local oscillator 
signals of the form 

V l 0 (0 = e i [“ lo(t) + 01l(t) ] 

(9) 

For Station #1 


v 2 0 (t) = e i [ u, 2o(‘) + 0i 2 (t)] 


( 10 ) 


Where co 10 andaj 20 are the local oscillator frequencies at the stations 1 and 2 
respectively. 0 u (t) and 0 12 (t) are the instrumental phases at the stations 1 
and 2 respectively. 


^20 ^ 


The mixed outputs at the two stations are of the form (assuming co 10 , 
v'(t) = Vj(t) V 10 (t) = e i[(^ - ^lo) 4 + 1 >T 1 (t) + 011 (t)] 


( 11 ) 


V 2 (t) = V 2 (t) V 20 (t) = e »[(“ - “20) t - + 0T2<O + ^2 (t) ] 

( 12 ) 


In one cross-correlation process, the signal from station #1 is offset by 
a model time-delay r d and the signals are multiplied together to produce a 
beat signal of the form 

V OUt= 

_ e i[(w - w 10 ) (t - T d ) + 0 X1 (t) + 0JJ (t)]. 

e i[(w - w 20 ) * " WT g + <t> T2(t) + (13) 
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where A</> x (t) = 0 x ,(t) - <*> x2 ( t) 

A0,(t) = 0„ (t) - 0 I2 (t) 

After the heterodyne operation, one can compute the correlation function which 
is proportional to the real part of the beat signal. 

R 0g, T d ) a £ Real (V out ) = A cos j^j (r g - r d ) + w, Q r d + A0 x (t) + A0,(t) ( 15 ) 

Let us assume that the two local oscillator frequencies are identical. Then 

R (7g,T d ) aC Ac0s C °( T g- r d) + % T d + A0 T (t) + A0,(t) ( 16 ) 

The correlation function has a maximum when r g = r d 
Then 

R max<>g’ T d) A C0S [ W o T g + A< M l ) + (17) 

If more than one local oscillator is used, co Q is the "effective" local 
oscillator frequency. 

Also because of the 2ttxv ambiguity involved in the inversion of trigono- 
metric functions, only the time-dependence of phase may be obtained from the 
correlation function. That is, the phase may be determined except for an 
unknown additive constant. The fourier- analysis technique is used to determine 
the fringe frequency. 

"c = w (* g -*d) + w o *d +a ^t + A ^i (18) 

y 

where the "dot" represents the time derivative of each parameter. 

Position determination depends directly on the measurement of time- 
delay, T g . The accuracy of r g in turn depends on having a broad -bandwidth 
source and a capacity for recording and processing over a sufficiently wide 
bandwidth in the cross correlation of signals. The significance of increased 
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bandwidth is reflected in the nature of the correlation function. For a mono- 
chromatic signal 

R(r ) » Acosfco 0 r d + w(r -r d ) + A 0 T + A<M 

(19) 

For a wide band signal, the cross-correlation function becomes 


R (V T d) = 


sin A co (T d - T g ) 

Aco ( r d 1 T g) 


cos[co 0 r g + A 0 T + A0j] 


( 20 ) 


Where Aw is the recorder bandwidth. The uncertainty in determining 
r g decreases with increasing bandwidth as (Aco)" 1 . Most recorders used in 
VLBI have a recording range from 360 kHz to 2.0 MHz. For bandwidths larger 
than 2.0 MHz, the bandwidth synthesis technique is utilized to cover a large 
bandwidth. In the bandwidth synthesis technique, the local oscillator and paramp, 
at both the terminals are switched in frequency in steps and synchronously 
retuned by means of programmed bias voltages. Careful phase calibration is 
done at the beginning of each switching. 

Why is it necessary to switch the received frequency over a wideband 
rather than to simply record a band sufficiently broad to yield high accuracy in 
differential delay? The basic reason arises from the lack of real-time wideband 
data link between the receiver elements of a VLBI system. The necessity to 
record the predetected signal causes.the instantaneous bandwidth of the VLBI 
system to be limited by that of the recorder, which in practice is a very small 
fraction of the received signal frequency. This constraint is overcome by 
varying the received frequency synchronously and coherently in cycles consisting 
of several steps. The bandwidth spanned is then not limited by the recording 
system but only by the bandwidth of the radio frequency amplifiers used. 

Thus, the following steps are undertaken in the determination of the 
fringe phase. 

1) RF signals are received, amplified and converted to video at as 
large a bandwidth as possible. 

2) At each terminal the video signal is clipped so that all positive 

areas of its signal have a value of 1 and all negative areas are zero. The 
clipped signals are sampled at the Nyquist rate. of twice the bandwidth and the 
bit strings are recorded on video or magnetic tape. Timing information is 
inserted in the data, format. . . ■ 

3) The tapes are brought together at a correlator which may be a digital 
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computer (CDC-6400, IBM-360-91, 75, 95) or a computer of special design. For 
each record, the bit streams are shifted relative to one another by an amount 
that corresponds to the theoretically predicted geometrical delay plus any 
fixed instrumental delay. This serves to produce "white" fringes. The records 
are then correlated. Usually, it is advantageous to divide the correlation inter- 
val into a number of sub-intervals in each of which the earth’s rotation rate can 
be considered a constant. 

4) The correlation function for each sub-interval or "frame" is 
corrected for sampling and for finite record length. Then it is multiplied by 
the trial function e j 5 1 where co T is the sum of all frequencies that 
tend to rotate the fringes (eg: earth's rotation rate, fringe-rotation rate and 
other unknowns) and 5t is the length of the sub-interval. This is repeated 
for a number of guesses for the unknown phase rate. 

5) Each of the modified correlation functions computed above is sub- 
jected to Fourier Transform. The case for which the transform peaks is taken 
to be the solution. For this case, the total phase-delay that has been applied to 
the signal is recorded. The wider the bandwidth, the fewer are the deep fringes 
and therefore the more precisely the white fringe region can be identified. 

For infinite record length, the correlation function has the form 

oo 

R (r) = J x(t) • x(t + r) dt (21) 

But in practice, for finite record length 

T + At/2 

R 00 to = / x(t - At/2) x (t + At/2) dt 

T - At/2 ( 22 > 

The purpose of applying a correction is to create a modified correlation 
function R Q (r) whose transform exists and is equal to the average value of the 
true transform. It is created by using a function known as the Hanning weighting 
function. It is also called the "Lag-window. " The function is characterized as 
follows : 


oo(0) = normality 

oj(-t) = oj(t) symmetry 

oo(r) >T =0 where T is the length of the record. 
n n 
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Two widely used functions are 



+ 



= 0 


;M<T n 


; M > T n 


( 23 ). - 


and 


oj ( t ) = 0.54 + 0.46 cos y 
= 0 


7 1 T 
m 


; r <T 


r l> T 


(24) 


It is usual to guess a number of trial-delay offsets and to evaluate the 
Fourier Spectrum of the fringes corresponding to each offset. The offset that 
causes the spectral response to peak is chosen to be the correct value of delay, 
i„e. the value that establishes the white fringe condition. 

Let us suppose that the signal received in one channel is delayed so that 
at that moment all frequencies within the bandwidth are in phase. If the ratio 
of BW to the nominal frequency is small, then all the frequencies will remain 
almost perfectly in phase during some subsequent time interval. Within this 
interval, the interferometer response will remain almost perfectly sinusoidal 
in shape, even though the bandwidth is of nonvanishing extent. Then the inter- 
ferometer is said to be operating in the "white fringe" condition, i.e. all the 
frequencies are present in the deep fringes. 

The power spectrum is denoted by 

T 

S(w) = limJ_ f R(r)e JWT dr 

T-> °o T j (^ 5 ) 

This in general will have the complex form 

S(w) = A(w) (26) 

Where A <p is the phase-delay evaluated at the beginning of the recording interval 
and gj t represents the total fringe rotation rate during that interval due to all 
causes. The interval is chosen small enough that is taken as a constant. 
Some effects that contribute to the to T are the fringe rate offset, rotational 
motion of the earth, and accumulated errors due to clock instability etc. If a 
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modified visibility function S' (<u) is defined as 

s'(co) = Y j* S(co)e j< V dt 


( 27 ) 


Where co' T is the trial function which includes a guess for the unknown term, 
then it is clear that S’ (to) is maximized when o/ t =oj> t ; because S (co) is close 
to a sinusoid and the integrated product of two sinusoids is nonvanishing when 
and only when their frequencies are equal. In practice, the offset can be deter- 
mined by this method to better than 0.001 Hz. Usually, the integrations with 
respect to r and t are more conveniently reversed, for then the fourier 
transform need not be evaluated so frequently. 


As mentioned earlier, maximum correlation is obtained if the recordings 
were offset by an amount equal to the time required for the wave to traverse the 
difference in distance between the source of the two sites. As the offset is 
changed from this amount, the two recordings would become more and more out 
of phase (up to a point) and the cross correlation would decrease. 


But the cross-correlation will be a maximum for many values of r 
(depending on the length of the recording) each of which differs from the other 
by the equivalent of 2tt in the phase of the signal. Thus 


cr = nX + SX ; 0 ^ SX < X 


(28) 


27rfr g = 27m + 0 

fr =„ + t <29) 

g 277 

Where n is an integer which is not determinable from a single frequency mea- 
surement. If we record a bandwidth /if instead, the ambiguities will now 
spread only over a time offset interval approximately equal to (Af)' 1 . Then, 
the cross-correlation as a function of the difference, At, between the offset 
and the time delay is proportional to 


sin 


2itMAt 


gj 


277 A f A t 


g 


COS (277fAT g ) 


(30) 


The first term provides the Sl " — envelope that restricts the ambiguity 
in the correlation maximum of cos 27ifAr g to lie within a width of ( Af ) -1 

There are alternate ways to minimize the error in the determination 
°f Ar g . The fringe phase is given by the expression 
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2nfr = 2 irn + <p 
g y 


( 31 ) 


then 


dip 

df 


= 2 nr 

g 


) dip 
T s 2v df 


(32) 


i.e., the delay sought is simply the slope of the phase versus frequency curve: 
at a given instant, the cross -correlation will have a maximum at a different 
phase for each different frequency in the sampled band. The time offset stays 
constant at all frequencies. In practice, a geometrical offset is used before 
correlation is done. Then the error in the estimation of r g is 

At = A<t> = A< ft 

& 2n(f - f . ) 2ttA[ 

v max rmiv ^ ( o o ) 

Where Af is the effective wideband (360 kHz, 2.0 MHz, 30 MHz, etc.) and A <p 
is the net phase change in the bandwidth interval. 7 - + A r g is the true time- 
delay at any instant. One can measure the phase/frequency slope to greater 
accuracy than is possible in extracting the bit offset from the correlation spec- 
trum where the region of ambiguity has a width of ( Af) -1 . 
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4. VL i r SYSTEM RESOLUTION 

The sources for the interferometric measurements can be either 
natural or artificial radio sources. Among the natural radio sources, the most 
desirable are the strong, broad-band point-sources of extremely small angular 
dimensions (a few arc seconds) and with negligible proper motion. These are 
quasars. 


Detectability of the noise signals from the tracked sources is a critical 
problem in interferometry. The total background noise coupled with the receiver 
noise is the dominant contributor to the system noise temperature. The effec- 
tive signal noise ratio of an interferometer pair is given by: 


S 

N 



T a 

Y 1 (BW) 7 
S 2 


% 


(34) 


Where T a and T a are the antenna temperatures of stations 1 and 2 
respectively, T s ^ and T s ^ are their respective system temperatures, (BW) is 
the instantaneous recorded bandwidth and t is the integration time. 


S o A l*l 
2k 

S 0 A 2 v 2 (35) 

2k 

Where S Q is the observed source flux density in Watts/m 2 cps, 

Aj and A 2 are the aperture areas of antennas 1 and 2 respectively, 
and t} 2 are their efficiencies, k is the Boltzmann Constant -- 1.38 x 10~ 23 
joules°K. 




f (AjfJj) (A 2 t? 2 ) (BW) (t) 



(36) 
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Where D ] and D 2 are the diameters of antennas 1 and 2 respectively. 

A few typical examples of signal to noise ratio is illustrated below: 

Rosman-Mojave Baseline t ATS-3 Satellite ] : 

Dj = 85 x 0.3 meters 

D 2 = 40 x 0.3 meters 

= V 2 = 0.5 

(BW) = 360 kHz 

r = 1 sec 

T = T = 65°K 

S 1 *2 

S « 2.77 X 10- 18 watts/m 2 The transponder output 
° is 446 watts. 

The signal to noise ratio will then be (For the ATS-3 Satellite Case): 

S/N “ 75 db 

Agassiz-Owens Valley Baseline [ ATS-3 Satellite ] : 

Dj = 84 x 0.3 meters 

D 2 = 130 x 0.3 meters 
V l = r? 2 = 0.5 
(BW) = 360 kHz 
t =100 sec 
T = 350° K 

;■ s t • 

T = 180°K 
S 2 

S Q = 2.77 x 10‘ 18 watts/m 2 
S/N' 85 db 

Similar signal to noise ratio computations can be made for the case 
when the observed sources are quasars. Let us assume that a quasar of flux 
density 10 units is monitored at the two terminals. 
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ThenS 0 = 10 flux units = 10 ' 25 watts/m 2 . 

For Rosman-Mojave baseline (integration time r s 100 sec) 

— ^ 16 db 

N 

For Agassiz-Owens Valley vaseline (integration time r~100 sec) 

| a 15 db. 

N 

The use of the satellite transponder as a radio frequency source permits 
high resolution time-delay interferometry because the spectral energy density 
of the transmitted signal is extremely high compared to celestial sources.; Flux 
densities correspond to an apparent source brightness about more than 100 db 
greater than the brightest quasar source at C-band. 

Time-Delay Resolution Uncertainties: Significant timing errors are due to 
three principal causes: 

!• clock epoch-alignment errors between the two stations.' 

2. long-term phase drifts in the atomic frequency standards used at 
the two stations. 

3. recording-bandwidth limitations. 

The clock-epoch alignment error (or the clock offset error) can be 
treated as an unknown while observing quasars and the time -delay observations 
from several sources will yield the clock offset parameter. While observing 
satellites, one has to resort to independent means for synchronizing clocks or 
accurately determining their relative offset. In our experiments, both the pro- 
cedures were adapted. Clock synchronization experiments were conducted 
using the dual transponders of the ATS-3 and ATS-1 satellites and the clocks 
were synchronized to within 10 nanoseconds. Also, extensive quasars observa- 
tions from the same sites also yielded the clock offset as a solution from the 
time-delay equations. 

Over sufficiently long observing periods, phase drifts in the atomic 
frequency standards also produce significant timing errors. Both rubidium- " 
vapor and hydrogen maser frequency standards are employed in the VLBI 
experiments. Recent data indicate that a stability figure of 1 part in lO' 12 is 
achieved with rubidium standards; hydrogen masers are about 1.5 to 2.0 orders 
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of magnitude better stability. Over a 24 hour period, the maximum accumulated 
time-error with a rubidium standard can be expected to be around 0.1 m sec. 


The time-delay 
the formula: 


resolution of the VLBI recording system is given by 


1 1 
6at “ 2tt(BW) /S/N ’ v 


(38) 


Where S/N is in the signal-to-noise ratio, (BW) is the recording badnwidth 
( ~ 360 kHz) and v is the number of independent observations. Thus the fac- 
tors affecting the precision of the VLBI measurements are the overall signal- 
to-noise ratio, the instantaneous recording bandwidth, the stability of the 
reference local oscillator and the unknown atmospheric and instrumental phase 
fluctuations. 


Based on the equation for e Ar , it is possible to compute the time-delay 
resolution for the case of satellites and quasars. A nominal recording band- 
width of 360 kHz is assumed for both the cases. 


For the ATS -3- Satellite Case: 
S/N » 75 db « 3 x 10 7 
v = 1 


Then, e A 0.1 nsec 

1 At 

In practice, however, the time-delay resolution that is achieved is 
limited by the frequency stability of the reference oscillator. The attainable 
time-delay resolution is of the order of 5 nanoseconds. 

For quasars ( ~ 10 flux units): 


S/N ^ 15 db 
v = 1 

e. ~ 75 nsec. 

At 

It is therefore possible to develop tne overall system resolution criteria 
for the Rosman -Mojave and Agassiz-Owens Valley baselines. The following 
table (Table 3) summaries the criteria. 
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TABLE 3 


VLB! System Resolution Criteria 


Criterion 

Parameter 

Satellite 

(ATS-3) 

Quasar 
(10 F.U .) 

Signal-to- 

Time-Delay 

5 nsec 

75 nsec 

Noise Ratio 

Fringe-Rate 

500-1000 fxHz 

5-10 mHz 

Stability of 
the Reference 

Time-Delay 

5 nsec 

10 nsec 

Frequency 
Standard „ 

- 1 Part in 10 

Fringe-Rate 

200 /x Hz 

1 mHz 


Expected System Resolution 


Source 

_ ATS-3 Satellite Quasar 


Parameter 

Time-Delay 

Fringe-Rate 


5 nsec 

100 - 500 /jtfz 


75 nsec 
5 mHz 
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5. DATA ANALYSIS 

5.1 OBSERVATION ROUTINE 

A basic observation routine was followed in which one station (usually 
Rosman) transmitted an uplink carrier signed with a preselected modulation. 
The recording bandwidth in the experiment was 360 KHz. Increased recording 
bandwidth improves the measured time -delay resolution. However, marked 
improvement can be achieved with narrowband (360 KHz) recording if signal 
of specific spectral distributions are transmitted/received to/from the space- 
craft transponder. Several modulation schemes were utilized as VLBI beacon 
signals to test their effect on the cross-correlation function. The modulations 
utilized in the experiment are listed in the following table. 


TABLE 4 


Signal Modulations 


Modulation 

Description 

Noise Translation 

Frequency translated broadband 
noise. Noise bandwidth 12 MHz 

FM Noise 

Frequency modulated carrier using 
0-12 MHz bandwidth Marconi 
noise generator 

Spacecraft Noise 

Front-end noise emitted by the 
ATS-3 transponder at 4178 MHZ 

PRS 

Pseudorandom sequency modulation 
using clock-rates of 500 k bits/sec 
1000 k bits/sec, 2 M bits/sec, 4 
M bits/sec 

Range Tones 

Ensemble of range tones simultane- 
ously used to modulate the carrier. 
The tones are 160 Hz, 4KHz, 

20 KHz, 100 KHz and 500 KHz 


The basic observation routine called for Rosman to transmit an uplink 
signal at 6301 MHz with a preselected modulation. The signal is transponded 
back to the earth from one of the ATS-3 transponders (Transponder #2, 6301 
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MHz uplink/4178 MHz downlink). The downlink signal is received at Rosman 
and Mojave. Both stations recorded the data simultaneously for 9 minutes. A 
variety of modulation schemes were utilized in a daily routine. Immediately 
before and after each VLBI recording, SHF ranging (R and R) was conducted 
simultaneously from Rosman and Mojave using the dual transponders on the 
spacecraft. The ranging duration was normally three minutes. The ranging 
data is utilized to compute independently the spacecraft orbit for comparison 
with the orbital elements computed from the VLBI data. It also provides rough 
estimates of time-delay (r ) and fringe-rate ( v ( ) for the intervening VLBI 
recording period. It would 8 also have been - ’esirable to conduct a systematic 
test on the effect of varying SNR on the precision of Tg andy f (VLBI parameters). 
This however turned out to be not feasible since the spacecraft transponders are 
hard limited and there is little flexibility in varying the output power of the 
transponders. 

A typical experiment schedule is shown in Table 5. 

5.2 CLOCK-SYNCHRONIZATION 

The VLBI clock at Mojave was set by a portable cesium clock and was 
also regularly monitored by a microwave relay signal both derived from the 
neighboring JPL/Goldstone Echo-Site. The Rosman clock was checked by 
monitoring Loran-C transmissions from Cape Fear, North Carolina. 

The synchronization between the Rosman and Mojave clocks was 
achieved by the simultaneous exchange of microwave signal pulses of 10 jusec 
duration and a very sharp rise-time via the dual transponders of the ATS-3 or 
ATS-1 satellites. A precision of 10 nanoseconds was accomplished using this 
synchronization technique. More details of this time-synchronization experi- 
ment can be found in the report "VLBI clock synchronization tests conducted 
via the ATS Satellites" by J. Ramasastry et. al„ NASA/GSFC, X-553-71-514, 
November 1971. 

5.3 ATMOSPHERIC PARAMETERS 

Both the troposphere and the ionosphere are significant sources of error 
in the measurement of the time-delay and fringe-rate parameters. The iono- 
spheric error can be large at L-band (1600 MHz) and even at S-band (2200 MHz). 
At C-band (4000 MHz) which is the frequency that was used in our experiments, 
the ionospheric influence is extermely small. The atmospheric error becomes 
large at low elevation angles. Table 6 summarizes the ionospheric and atmos- 
pheric time-delay errors at various zenith angles. 


38 



TABLE 5 


First 

Session 


Typical Daily Experiment Schedule 


Stations: Rosman & Mojave 


Source: ATS-3 

Statellite 

Time 

Duration 

IVansition 

Stations 

Mode of Operation 

Comments 

(GMT) 

’ min ) 

[ min ) 




T + 0 

10 

5 

R/M 

Time synchronization test 

Dual SHF 

transponders 

used 

T+15 

3 

2 

R/M 

(R) - VHF Ranging \ 
(M) - SHF Ranging / 

Two frequency 
ionospheric 

T+20 

3 

2 

M/R 

(M) - VHF Ranging 1 
(R) - SHF Ranging / 

calibration 

T+25 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T+30 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
noise trans- 
lation 

T +40 

3 

2 

R/M 

R/M 

SHF Ranging 
VLBI Recording 

Simultaneous 
Rosman Xmit. 

T +45 

6 

4 

noise trans- 
lation 



T +55 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T+60 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
FM noise 

T + 70 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T + 75 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
FM noise 

T + 85 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T+90 

6 

4 

R/M 

VLBI Recording 

No Xmission 
Spacecraft 






noise 

T + 100 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T+105 

6 

4 

R/M 

VLBI Recording 

No Xmission 

Spacecraft 

noise 

T + 115 

3 

2 

R/M 

SHF Ranging 

Simultaneous 
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TABLE 5 (cont.) 



Time 

(GMT) 

Duration 
( min ) 

Transition 
( min ) 

Stations 

Mode of Operation 

Comments 


T + 120 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit 


T + 130 

3 

2 

R/M 

SHF Ranging 

PRS 500 kbit/ 
sec 

Simultaneous 

First 

T + 135 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 

Session 

T + 145 

3 

2 

R/M 

SHF Ranging 

Prs . 500 kbit/ 
sec 

Simultaneous 


T + 150 
T + 210 

10 

60 

5 

R/M 

Intermission 

Time synchronization test 

ATS- 3 dual 


T + 225 

3 

2 

R/M 

SHF Ranging 

transponders 

used 

Simultaneous 


T +230 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 


T + 240 

3 

2 

R/M 

SHF Ranging 

noise trans- 
lation 

Simultaneous 


T +245 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 

Second 

Session 

T + 255 

3 

2 

R/M 

SHF Ranging 

Noise trans- 
lation 

Simultaneous 


T +260 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 


T + 270 

3 

2 

R/M 

SHF Ranging 

Noise trans- 
lation 

Simultaneous 


T + 275 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 


T + 285 

3 

2 

R/M 

SHF Ranging 

FM noise 
Simultaneous 


T + 290 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 


T + 300 

3 

2 

R/M 

SHF Ranging 

FM Noise 
Simultaneous 


40 












TABLE 5 (cont. ) 


Time 


Transition 

Stations 

Mode of Operation 

Comments 

, (GMT) 


(min)-,', . 



T +305 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
FM noise 

T + 315 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T + 320 

6 

4 

R/M 

VLBI Recording 

No Xmission . 




Spacecraft 



- 



noise 

T + 330 

3 

2 

R/M . 

SHF Ranging 

Simultaneous 

T +335 

6 

4 

R/M 

VLBI Recording 

No Xmission. 
Spacecraft 






noise 

T + 345 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T + 350 

6 

4 

R/M 

VLBI Recording 

No Xmission. 
Spacecraft 






noise 

T + 360 


60 


Intermission 


T +420 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T +425 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
Noise trans- 
lation 

T +435 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T +440 

6 

■ 4 

R/M 

VLBI Recording 

Rosman Xmit. 
Noise trans- 
lation 

T +450 

3 

2 

R/M 

SHF Ranging 

Simultaneous 

T +455 

6 

4 

R/M 

VLBI Recording 

Rosman Xmit. 
Noise Trans- 
lation 


T +465 

3 

2 

R/M. 

SHF Ranging 

Simultaneous 

T +480 

3 

2 

R/M 

(R) - VHF Ranging \ 

(M) - SHF Ranging / 

Two Frequenc 

T + 485 

3 

12 

R/M 

(r) SHF Ranging ( 

Ionospheric 




R/M 

(M) - VHF Ranging / 

Cal ibration 

T + 500 

15 

— 

Time Synchronization Test 

Dual trans- 

T + 515 

•1 ■ 


1 

End 

ponders on 
ATS-3 used 


Y 









TABLE 6 


Typical Values of Refractive Bias 

(Ref: Refractive Corrections in High- Accuracy Radio Interferometry: D.F. 
Dickinson, M.D.Grossi &M. R. Pearlman, J. Geophy. Res. 75 (8), 1619-1622 
March 1970) 


Zenith Angle 

Medium 

X-band 
10 GHZ 
(meters) 

C-band 
5 GHZ 
(meters) 

L-band 
1.5 GHZ 
(meters) 

0° 

Atmospheric 

-2.3 

-2.3 

-2.3 


Ionospheric 

0.1 

0.3 

3.3 


Total 

-2.2 

-2.0 

1 .0 

45° 

Atmospheric 

-3.2 

-3.2 

-3.2 


Ionospheric 

0.1 

0.4 

4.4 


Total 

-3.1 

-2.8 

1 .2 

90° 

Atmospheric 

-86.6 

-86.6 

-86.6 

(Horizon) 

Ionospheric 

0.2 

1 .0 

1 1 . 0 


Total 

-86.4 

-85.6 

-75.6 


Atmospheric data from various sources were collected for correction of 
propagation effects on VLBI observables. The temperature, pressure and 
humidity were measured at both stations with hydrothermograph and microbaro- 
graph and monitored continuously by strip-chart recordings. These data are 
used to calculate surface refractivity and then to predict the tropospheric 
correction along the line of sight. 

Several means were employed to obtain ionospheric data. A very 
accurate technique was utilized during our VLBI experiment. On several occa- 
sions during the VLBI experiment each day, the electron content along the 
line-of-sight was determined directly by conducting a two-frequency differential 
delay measurement via the ATS-3 satellite. This was accomplished by per- 
forming in rapid sequence SHF and VHF ranging from the two stations. Only 
a small extrapolation was needed to account for the non-simultaneity of the two 
ranging operations and bring both of them into a common time-frame. The 
difference between the VHF and SHF range measurements at any particular 
time gives the ionospheric propagation delay. This technique, though useful 
and precise was not always accessible. Furthermore, observations at Mojave 
were regularly made during the later afternoon when the electron content 
changes rapidly. It was therefore deemed desirable to have other ionospheric 
measurements. Faraday rotation measurements of ATS-1 and 3 VHF signals 




were acquired at geographically favorable locations by us and other experi- 
menters during the VLBI experiment. The Faraday polarization angle is 
directly related to the line of sight total electron content. 

Table 7 gives details about the tropospheric and ionospheric data 
collected during the experiment. 


TABLE 7 


Atmospheric Data Monitored During the VLBI Experiment 


Medium 

Parameter 

Method/ 

Instrumentation 

Site(s) 

Space- 

craft 

Troposphere 

Humidity 

Temperature 

Pressure 

Hydrothermograph 

Microbarograph 

Rosman 
Mojave , 


Ionosphere 

Actual Ionospheric 
Delay 

Two Frequency 
(VHF/SHF) 
Simultaneous 
Ranging 

Rosman 
Mojave ■ 

ATS-3 


Total Electron 
Content 

FHF Faraday 

Rotation/Polari- 

meter 

JPL/Goldstone 

Stanford 

Rosman 

Wallops- 

Island 

ATS-1 , 
ATS-3 
ATS-3 
ATS-3 

ATS-3 


Ionospheric 

Differential 

Delay 

Global Ionospheric 
Model ing/ Ground 
Based lonosonde Dat 
Dr. Rodney Bent 
DBA Systems Inc. 

3 



Figure 12a shows the time-delay corrections as applicable to a frequency 
of 4178 MHz for the Rosman-Mojave baseline for the quasar 3C-273 for June 3, 
1971. Shown are the ionospheric and tropospheric delay corrections. Also 
shown in the figure are the universal time versus elevation angle on the star 
for both Mojave and Rosman for June 3, 1971. Figure 12b gives the multiplica- 
tive factor for computing the ionospheric correction for any other frequency 
larger than 1 GHz. 
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TIME DELAY CORRECTION ( nanoseconds ) 


DATE= 3 JUNE 71 


FREQUENCY = 41 78 MHz 



Figure 12a. Ionospheric and Tropospheric Corrections for the Quasar 
3C-273B Observations from Rosman and Mojave 
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FREQUENCY GHz 
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2 . 


2 . 


2 . 


2 . 



Figure 12b. Multiplication Factor for Ionospheric Corrections 
for Frequencies > 1 GHz 

Figure 12c shows the time-delay corrections as applicable to a frequency 
of 4178 MHz for the Rosman- Mojave baseline for the ATS-3 satellite for May 26, 
1971. Also shown in the figure are the elevation angles of the satellite from 
both Mojave and Rosman. Figures 12a and c show sample results. The ionos- 
pheric and tropospheric time-delay corrections are stored in the computer for 
correcting the VLBI data. 

5.4 SHF RANGE & RANGE-RATE TRACKING 

An independent check on the accuracy of the VLBI parameters 
[namely, time-delay, r g , and fringe- rate, v { ] can be obtained from C-band 
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TIME DELAY CORRECTION ( nanoseconds) ELEVATION (deg) 


FREQUENCY = 4178 MHz 


DATE =26 MAY 71 



UNIVERSAL TIME (hrs) 


Figure 12c. Ionospheric and Tropospheric Corrections for ATS-3 
Observations from Rosman and Mojave 
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tracking data. The VLBI parameters have a simple relation with the range and 
rate measurements made at the same frequency. They are given by 

R - R 

^Mojave ~ Rosman 

T g = c (39) 

and A _ r 

^Mojave Rosman (401 

"f X 

where c is the speed of light and X is the wavelength of the radio signals. R 
and R denote the ground to spacecraft range and range-rate. Determining r g 
and v by R and R tracking is also of value in approximately predicting the 
VLBI 1 parameters. This saves computer time in the initial search process, 
even though it is not essential for processing VLBI data. SHF range and range- 
rate measurements are also used to compute the orbital elements independently 
for later comparison with the orbital elements derived from the VLBI data. It 
is for the various reasons cited above that R and R measurements from both 
Mojave and Rosman were incorporated in the observation routine (Table 5). 

5.5 VLBI TAPE PROCESSING PROCEDURES 

The recording of the high bit-rate (~ 720 kbits/sec) is done by using 
the Ampex VR-660C helical scan video tape recorders which have a high density 
bit packing capacity. The reels of video tape used in the experiment store up 
to 96 minutes of observational data (~720 kbits/sec) as against individual 2400 
feet 9-track tapes which can only record 3 minutes worth of the same data. 

At the conclusion of the observation phase of the experiment, the tapes 
along with the recorders were transported to the Central processing facility 
where the playback system is used to transcribe the data on video tapes on to 
9 track digital tapes. The cross -correlation and fringe -fitting is accomplished 
in the general purpose high-speed computers, the IBM 360/95 at the Goddard 
Space Flight Center and the CDC-6400 at the Smithsonian Astrophysical 
Observatory. 

Seven computer programs are available to process the data. First is a 
pre-processor which does a quality check and realignment of the nine-track 
data tapes and transcribes the data to seven -track tapes. One of two fiinge 
processors is used to do the correlations and determine delay and doppler 
for each pair of tapes. One fringe processor is used on satellite data, the other 
on stellar data. Also in use are two pairs of least-squares analysis programs. 
Each pair consists of a program which uses delay data and the other one which uses 
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delay-rate data. One pair is for stellar data and the other is for satellite data. 
Brief descriptions of all these programs are provided below. More details 
about these programs may be found in the publication "Documentation of Compu- 
ter Programs Used in Very Long Baseline Interferometry Experiments," 

S. Ross and Company, November 1971, NASA Contract NAS-5-20247. 

5.5.1 Garbage Collection Program - Tapes 


The digitized data is recorded at the observing stations on Ampex VR-660C 
helical-scan videotape recorders. These video tapes are played back and the 
data is transferred to nine-track computer tapes. Along with the data on both 
video and computer tapes are recorded timing and syncrhonization codes. The 
data on the nine-track tapes are processed to align these codes on the proper 
word boundaries, inserting or deleting bits where necessary. This operation is 
performed by the Garbage Collection program. Incidental to this operation is 
the transcription to seven-track tapes. This is done because the SAO computer 
cannot handle pairs of nine-track tapes. The nine track tapes are used in GSFC 
computers. 

5.5.2 Stellar Fringe Processor - VLBI1 


The Fringe Processor does a two dimensional search for a peak fringe 
amplitude. The inputs consist of parameters describing the searches in the 
time-delay offset domain and the fringe frequency offset domain. Also required 
is the integration interval and certain parameters related to the recording 
system. The brunt of the computing is spent doing the bit-by-bit cross-correla- 
tion of the two data strings. The program computes delay and delay rate from 
the VLBI geometry and shifts one of the data strings to compensate. As the 
integration proceeds, the predicted phase difference is adjusted by the predicted 
doppler frequency. These compensations cancel out the known offsets. The two 
dimensional search around these predicted offsets allows the analyst to ascer- 
tain errors in the predictions. These errors are used at a later stage in the 
computations. 

5.5.3 Satellite Fringe Processor - VLBI2 


This program operates in much the same way as VLBI1 except that the 
time delay due to the geometry is not computed by the program, but rather is 
input in microseconds. 
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5.5.4 Stellar Least Squares Analysis Program for Delay Data - POSTPR. 


This program does an analysis of the time delay errors from the Stellar 
Fringe Processor into corrections of source and baseline declination, baseline 
length, clock error, and baseline hour angle minus source hour angle. At least 
three sources must be used with at least three observations of each at different 
times of the day. The Least Squares technique is applied at two stages in the 
program, once on the data from each source and once on the results from the 
first stage. 

5.5.5 Stellar Least Squares-Analysis Program for Delay Rate Data - POSTDR. 


The Least Squares technique is used on the corrections to the predicted 
doppler frequency as ascertained by the stellar Fringe Processor. The assump- 
tion is made, for purposes of this program, that there are no errors in source 
coordinates. The data from at least three observations of natural sources is 
input and corrections to the hour angle and equatorial component of the baseline, 
as well as the clock drift rate are computed. An error propagation analysis has 
been added to this program to compute the likely errors in the corrections. 

5.5.6 Satellite Least Squares Analysis for Delay Data - DELAY. 


A sinusoid is fit to the delay offset data using the Least Squares technique. 
This program also produces a CALCOMP plot of the curve superimposed over 
the data points so that a visual check of the quality of the data can be made. Nine 
unknowns are estimated. 

5.5.7 Satellite Least Squares Analysis Program for Delay-Rate Data - DELRT. 


This program is analogous to the previous one except that the sinusoid 
is fit to the doppler offset data. Three unknowns are estimated. 

Figure 13 is a schematic diagram which outlines the data processing 
procedure. In processing satellite data, it was found sufficient to integrate only 
one record-length of data ( = 200 msec) to obtain unambiguous time-delay and 
fringe-rate values. However, the number of records integrated was increased 
to five (200 msec x 5 = 1 sec) since the sampling rate is identical to that of the 
range and range-rate data. The following table (Table 8) illustrates the fact 
that there was not much to gain by integrating more records. This is due to the 
fact that the signal-to-noise ratio (SNR) is so high in the case of satellites that 
one record of integration is sufficient to achieve strong correlation. 
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TABLE 8 


VLBI Parameters as a Function of Number of Records Integrated 
Satellite ATS-3 Run No. 159:22:11 (Year 1971) 


Records Integrated 
Over 

Doppler Search 
Interval (Hz) 

Fringe Frequency 
at Peak Correlation 

Maximum 

Correlation 

Amplitude 

1 

0.005 

7.855 

0.399981 

5 

0.005 

7.855 

0.398210 

10 

0.005 

7.855 

0.420594 

20 

0.005 

7.850 

0.411054 

50 

0.005 

7.850 

0.407456 
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Figure 13. VLBI Data Processing Scheme 


5.6 FRINGE PROCESSING RESULTS 


Four types of signal modulation were employed during the experiment 
session. The objective was to determine whether any particular modulation 
scheme would be especially advantageous for VLBI application. The modulation 
types were : 

1. Spacecraft Noise. Constant density analog noise transmitted from 
the spacecraft transponder front-end 

2. Noise Translation. Flat -noise spectrum having 12 MHz bandwidth 

3. Sidetone ensemble. A harmonic sidetone ensemble, using all of 
the ATS range tones turned on simultaneously. 

4. Pseudo-random Sequence. A pseudo -random sequence code using 
2047 sequences repeated at rates of 500 kbits/sec, 1 Mbits/sec and 
2 Mbits/sec. 
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Figures 14-22 show the results from processing sample data from typical 
runs of each type. The "a” figures show the correlation function amplitude and 
phase as a function of number of delay bits, where the center of abscissa scale 
(usually between points 20 and 21) indicates zero relative excess shift. Each "a" 
plot is for a specific doppler frequency. Each bit is of length 1.4 seconds. 

This bit offset is in addition to the geometric delay offset which corresponds to 
the zero offset (between points 20 and 21) in these figures. The additional 
offset bits where the correlation function maximizes are added algebraically to 
the geometric delay to get the true time-delay, r g . However, this technique 
gives time-delay resolution equivalent to the bit-length which is 1.4 ^ sec. An 
alternate technique which involves taking the slope of the phase versus frequency 
curve of the power spectrum gives a resolution for time-delay which is 2 orders 
of magnitude better than that obtained by using integral number of bits from the 
correlation plot. This technique will be discussed later on. 

The "b" figures show the amplitude and phase of the power spectrum 
corresponding to the "a" figures. The amplitude spectrum is characteristic of 
the power spectrum of the received signals and is a reliable indicator of true 
correlation. The phase versus frequency plot is the most important plot of the 
set. It gives indications of phase noise in the system (random phase excursions 
and jumps). 

The "c" figures show the "fringe spectrum." Shown in the figure are the 
peak correlation and power amplitudes as a function of the doppler frequency. 

The "a" and "b" figures are for one specific doppler frequency and each have a 
corresponding peak correlation and power amplitude. The "c" figure shows the 
distribution of the peak amplitudes with a doppler scan. One may see the -^1^- 
distribution in doppler space which is in line with the theory. The doppler 
frequency at which the peak correlation and power amplitudes maximize is the 
true doppler frequency (or fringe-rate) for that run. 

The slope of the power spectrum phase curve (i.e. "fringe phase") can 
be used to determine with great precision the excess shift required to produce 
a peak in the correlation function. It is done in the following manner. 


Total fringe phase, 


and. 

Therefore 


A0 = (A<P 0 +A<t> E ) = f 
A0 O = fAr 0 
A0 e = fAr E 


Ar 0 + At £ 


(41) 
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Where A<£ is the total fringe phase, A$ 0 is the geometrical fringe-phase, 
A is the excess fringe phase shown in the power spectrum phase plot, f 
is the RF frequency, A t q is the constant geometric -delay and At e is the 
excess time-delay that one needs to determine accurately. The total delay 
Ar = At q + At e is then easily computed. 


From the above expression: 


d (A0 e ) 
df 



(42) 


Where a denotes the slope of the fringe-phase curve. In other words, the 
slope of the phase curve is equal to the relative excess time-delay ("excess" 
here refers to the amount in variance with the geometric delay. The variance 
may be positive or negative.) This is important because one can measure the 
phase/frequency slope to greater accuracy (one to two orders of magnitude 
improved resolution) than is possible by extracting the bit-offset from the 
correlation amplitude curve. For instance, if the fringe-phase traces n cycles 
linearly across a frequency band of m Hertz, the relative excess-delay would 
be n/m seconds. A least squares fit gives this slope with a standard deviation 
( ~ resolution) of the order of 10 nanoseconds. The correlation amplitude 
curve will show a maximum around this computed excess-delay but at integral 
number of bits and the resolution is limited by the bit length which is 1.4 Msec 
for the present case. So, the correlation plots are never used to compute the 
excess delays. They are only used to check the sign and magnitude of the 
computed excess-delay. 


As is evident from the accompanying figures, each of the modulations 
used during the experiment exhibited its characteristic forms. Figures 14, 15, 

16 and 17 are examples of the noise translation and spacecraft noise modulations 
on the received signal. In the power spectra corresponding to these cases 
(Figures 14b, 15b, 16b and 17b), the shape of the bandpass filter and the flat 
noise spectra can be seen. The spectral phase curves are remarkably linear 
which is indicative of the extraordinary time-delay resolution. Figures 14 and 

16 are the cases where the bit offset is relatively large which is evidenced by 
the number of cycles in the corresponding spectral phase plots. Figures 15 and 

17 are for cases where the bit offset is not more than one bit. This is evidenced 
by less than a cycle of phase change in the spectral phase plots. Figures 15c 
and 17c show the typical Sin x profiles for the "fringe spectrum." 

x 

Figures 18, 19, and 20 are typical of the pseudorandom sequence modu- 
lation. Figures 18a and 18b are the correlation function and power spectrum 
respectively of the 500 kbit/sec rate PRS modulation. The sixbit width of the 


53 





J 


159-22-11 ATS 3 N-* 

SEARCH FREQUENCY NUMBER i 


SEARCH FREQUENCY 

7.9500 

maximum AMPLITUDE 

.428865 


NO, 

AMP 

phase 

amplitude 





RHASE 







-180 

-90 

0 90 

180 

l 

.003922 

172.693 



I 

I 



*1 

2 

.006971 

“153.251 



I 

1 * 



T 

3 

.020128 

1 1 5 . 099 

I* 


T 

T 


0 

I 

4 

.009819 

’ .422 

l» 


I 

I 


0 

I 

s 

.019031 

1-2.220 

I* 


I 

I 


0 

I 

6 

.032064 

19.442 

1 * 


I 

I 


0 

I 

1 

.010132 

“148.264 

I* 


I 

I * 



I 

o 

,058240 

27.234 

1 * 


t 

I 


0 

I 

9 

.062953 

-67.807 

1 • 


I 

1 


0 

I 

n 

.039215 

41.70* 

1 « 


I 

I 


0 

I 

1 1 

.155147 

-67.358 

I * 


I 

I 


0 

I 

12 

.0«3766 

-158.128 

I * 


I 

I * 



I 

13 

• 240 349 

-7". 620 

I * 


I 

I 


0 

I 

14 

.690930 

-165.850 

1 

0 

I 

I* 



I 

IS 

1.000000 

102.913 

I 


• I 

I 


0 

I 

IS 

.722-45 

7.271 

1 

0 

I 

I 


0 

I 

IT 

.282096 

-83.905 

I « 


I 

I 


0 

I 

i e 

.063639 

-16.908 

1 « 


I 

I 


0 

I 

19 

.1*0250 

-101.113 

I * 


I 

I 

« 


I 

?•> 

.0*23 .7 

16. .353 

I « 


t 

I 



<* I 

21 

,058419 

-95.933 

I * 


t 

I 

« 


I 

2? 

.054772 

169.380 

I * 


I 

I 



* I 

23 

.011940 

-48.049 

I* 


I 

I 


0 

I 

24 

.035542 

158.253 

I * 


I 

I 



* I 

2S 

.010442 

54.846 

I* 


l 

I 


0 

I 

26 

.014000 

-154.836 

I« 


I 

I • 



I 

27 

.015611 

66.905 

I* 


I 

I 


0 

T 

2 8 

.001568 

-161 ,210 

ft 


I 

I * 



I 

29 

.010474 

69.474 

I* 


I 

I 


0 

I 

3' 

,003992 

-67.961 

» 


I 

I 


0 

I 

31 

.006676 

151 .284 

« 


I 

I 



• I 

3? 

,008922 

-49.933 

I* 


I 

I 


0 

I 

33 

,002621 

-145.616 



I 

I * 



I 

34 

.003894 

75.828 

* 


I 

I 


0 

I 

35 

.003184 

-8.214 

» 


I 

I 


0 

I 

36 

.002928 

1 1 .870 

» 


I 

I 


0 

I 

37 

.003464 

61 .397 

* 


I 

I 


0 

t 

38 

.002961 

102.351 

« 


I 

I 


0 

I 

39 

.000525 

-97.606 

4 


I 

I 

« 


I 

4,1 

,0027-9 

-1 -5.427 

« 


I 

I 

0 


I 


Figure 14a. Correlation Function Versus Delay Offset for Noise Modulation - 
Case of Large Delay Residual 
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Figure 15c. Fringe Spectrum for Noise Modulation 
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Figure 15c. Continued 



159-22-27 S/C NO i SF 


SEAoCH fPEOUENC T NiiMHtR IS 
SEAoCH EPEOUENC* 6-7U0O 

MAXIMUM AMPLITUDE .378901 


NO. 

AMP 

PHASE 

AMPL ITUCIE 





PHASE 









-180 

-90 0 

90 

180 







. 




----* 

1 

. '05928 

-87.367 

« 



I 

I 

0 


I 

2 

.J13ARA 

-144.798 

i« 



1 

t • 



T 

3 

. A027ft7 

lpS.553 

a 



I 

I 


0 

I 

A 

.• 25837 

-90.030 

I « 



I 

I 

0 


I 

5 

. 31533 

188.355 

I * 



1 

I 



* I 

6 

.014129 

-98.433 

t* 



1 

1 

0 


I 

7 

.075021 

172.915 

I » 



I 

I 



*1 

8 

. "48968 

72.52: 

T • 



I 

I 


0 

I 

9 

. 75683 

174.799 

I « 



I 

I 



•I 

10 

.171214 

75.522 

I 4 



1 

I 


0 

T 

11 

.,'150871 

- 1 7.602 

I * 



1 

T 

0 


I 

12 

.33A176 

71.625 

I « 



I 

I 


0 

I 

13 

.814730 

- 21 .660 

I 


0 

I 

I 

0 


I 

1 A 

l.wOOOOU 

-116.873 

1 



* I 

I 

0 


I 

15 

, 684444 

1 52 • 1 'U >> 

1 

0 


I 

I 



» I 

1 ft 

.211986 

6 O .201 

I 



I 

I 


0 

f 

17 

,.,9801 1 

136.686 

I 4 



I 

I 



* I 

18 

.1367 2 

43.184 

r • 



I 

I 


0 

1 

19 

.019328 

-68.928 

1 * 



I 

1 

0 


I 

20 

.,'70235 

46.986 

1 * 



I 

I 


0 

I 

21 

. , 43460 

-59.917 

I « 



I 

I 

0 


I 

22 

.19511 

61.633 

1 * 



I 

I 


0 

I 

23 

. 33386 

- 46 .667 

I 4 



I 

I 

0 


I 

2A 

.11199 

140.805 

T* 



I 

I 



* I 

25 

1879 1 

-30.934 

I« 



I 

I 

0 


I 

26 

.17699 

-124.836 

I* 



I 

I 

0 


I 

27 

. '0 3661 

177.727 

* 



I 

I 



0 

28 

. 02653 

-146.989 

0 



I 

I * 



I 

29 

.,.04883 

84 ,6p5 

« 



I 

I 


0 

I 

30 

."13853 

- 1 n4 . 1 82 

I* 



I 

I 

0 


I 

31 

.003468 

1 1 3 . 46 1 

0 



I 

I 


0 

T 

32 

.•033'. 3 

-24.857 

« 



I 

I 

0 


I 

33 

.006086 

179.533 

* 



I 

I 



* 

3A 

.003536 

13.777 

0 



I 

1 

0 


I 

35 

.003886 

83.109 

0 



I 

I 


0 

I 

36 

• 1 06995 

111. 881 

0 



I 

I 


0 

I 

37 

.002279 

-24.710 

0 



I 

I 

0 


I 

38 

.001884 

-17.112 

0 



I 

I 

0 


I 

39 

.<-'02234 

123.249 

0 



I 

I 


0 

1 

AO 

.005662 

-179.713 

0 



I 

« 



I 


Figure 16a. Correlation Function for Broadband Spacecraft Noise - Case of Large Delay Residual 
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Figure 17a. Correlation Function for Broadband Spacecraft Noise - Case of Small Delay Residual 
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Figure 17c. Fringe Spectrum for Broadband Spacecraft Noise 
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Figure 18a. Correlation Function for Pseudo Random Sequence Modulation (500 kilobits/sec) 
Case of Large Delay Residual 
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Figure 18b. Power Spectrum for Pseudo Random Sequence Modulation (500 kilobits/sec) - 
Case of Large Delay Residual 
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Figure 19a. Correlation Function for Pseudo Random Sequence Modulation (1 Mbit/sec) - 
Case of Large Delay Residual 
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Figure 19b. Power Spectrum for Pseudo Random Sequence Modulation (1 Mbit/sec) 
Case of Large Delay Residual 



159-21-23 ATS-3 P-S 


SEARCH FREQUENCY NUMdER 
SEARCH FREQUENCY 
MAXIMUM AMPLITUDE 


19 

It. 0TC909Q 
0 .34 1 8C543935064 


<1 

o 


► 


AMPL I TUOE 


+ ♦ 


-180 -90 

- ♦ ♦ — 


PHASE 

0 


90 ISO 


\ 

9 • 9 053 3 r 

21 .742 

* 

2 

P.022966 

-173.527 

| * 

3 

C. 026937 

44 . 163 

I 

A 

9.01 91 36 

-119.170 

I * 

5 

0.022o6E 

121 .727 

I * 

6 

0, 921 59f 

-40 .654 

I 4 

7 

0.918967 

157 . 166 

1 * 

8 

O. P3293F 

1 1 .468 


9 

9 . ° 331 86 

-1 34 .941 

I 

1 9 

9. C 15732 

68. 975 

I * 

1 1 

3.01682C 

-130.922 

I * 

1 2 

0.057856 

72.739 

I 

1 3 

0.12 02 99 

-79.060 

I 

1 A 

0. 1 43576 

124 .793 

I 

15 

0. 16 9A4Z 

-45.689 

I 

i e 

0. 1591 40 

161. 064 

I 

l 7 

0 . 27592F 

1 .753 

I 

1 8 

9 • 6 A 37 72 

- 1* 2. 45? 

I 

1 9 

0.991686 

87 . 385 

I 

20 

1 . n 0 r. n 0 o 

-36.297 

I 

21 

0. 71 0515 

-156.044 

I 

22 

0. 436762 

69.50? 

I 

23 

9.251915 

-91 .635 

I 

24 

O • 1 8625* 

105.364 

I 

25 

0 . 1 3678! 

-64 . 323 

I 

26 

0. 1 009^8 

141. 266 

I 

27 

C. 068373 

-4.411 

I 

28 

O .043457 

-161 .629 

I 

29 

0.02794f 

—6 .477 


30 

9 . C‘ 3 A 7 87 

1 75. 807 

I 

31 

0.038999 

46 .378 

I 

32 

0 . 027640 

-9e. 928 

l 

33 

0.021176 

92.525 

I * 

34 

0. 022326 

—66 .993 

I 4 

35 

0.01 6716 

168. 189 

I 4 

36 

0.00 7504 

57. 029 

* 

37 

0 . 0 1 2940 

-11 3.737 

I * 

38 

0.00o559 

82.197 

* 

39 

0.008551 

53 .960 

I 4 

40 

0.014A83 

-1 1 3.661 

I 4 


Figure 20a. 


Case of Small Delay Residual 



PLOT OF FREQUENCY SPECTRUM • PEAK VALUE 


1 .655396 


NO. 

AMP 

PHASE 

A MPL ITUDE 

rMAo 

-180 -90 0 




+ 

II * 

1 

0.007^51 

-24.644 

* 

II 4 

2 

0 • 0 09.352 

-21 .905 

I* 

11 4 

3 

0.011 337 

-26.469 

14 

I I 4 

A 

0. 01 26*e 

-28.416 

I * 

II * 

5 

0. 01 9968 

-14 .061 

I* 

II * 

6 

0.035527 

-1 3.596 

I * 

1 I 4 

7 

0.053862 

-22. 199 

I ♦ 

II 4 

8 

0 .069934 

-31 • 375 

I * 

I I 4 

9 

0.084192 

-36.543 

1 m 

I I 4 

l 0 

0. 100547 

-40.632 

I * 

II .4 

1 1 

•). 1 2 3 65! 

-49.600 


I I 4 

1 2 

0. 1 54869 

-6 0.718 

I * 

II 4 

1 3 

0. 137013 

-67.827 

I 4 

II 4 

1 A 

0 .2135 75 

-71 . 51 1 

I 4 

I I 4 

1 5 

0.2321 55 

-76.222 

» 4 

II 4 

1 6 

0.2531 99 

-83.849 

I * 

II 4 

1 7 

0.26S*6f 

-91.232 

I * 

II 4 

1 8 

0. 331706 

-96 .467 

I 4 

II 4 

1 9 

0.391 605 

-10). .008 

I * 

II * 

? 0 

0.A59957 

-10 8.647 


II 4 

21 

0.517920 

-115.676 


II 4 

22 

0 • 554 A 84 

-121 .842 

I 4 

II 4 

23 

0.57458! 

-127.532 

i * . 

11 4 

24 

0. 595013 

-133.387 

I 4 

II * 

25 

0.61 9941 

-139.306 

r * 

II * 

26 

0 .64451 0 

-1 45. 065 

I * 

II 4 

27 

0.660087 

-150. 692 

| * 

T I * 

2 5 

0 .746727 

-155.976 

I 4 

# 

I I 4 

29 

0.375679 

-160.936 


A I I 4 

30 

0. 9051 09 

-166.261 

I 

* [ J * 

31 

932374 

-172.437 

I 

* 14 

3? 

0 • 944505 

-178 . 951 

I 


33 

0 .974.503 

175.176 

I 


34 

l . ooooor 

169.613 

I 

4 II 

35 

0.956482 

163 . 40 1 

I 

I t 

3 6 

0. 80 3746 

156.799 


t i 

r | 

37 

0.595240 

151 .484 

j * 

J J 

38 

0.384950 

1 49. 1 76 

i * 


39 

0 . 228292 

147.179 

I 4 

I I 
| | 

40 

0 .*1665 11 

136.326 

-I * 

X X 


LEAST SQUARES FOR PHAS- (IN CYCLES) VS FREQUENCY (IN HZ) 

PHASE = 1902807690- )5) * NO. 0. 46*3599540-0! 

SIGVA IN PHASE = 9.7213377390-02 
SIGMA IN SLOPE = 9. 1 53*5? *980-07 

Figure 20b. Power Spectrum for Pseudo Random Sequence Modulation (500 kilobits/sec) 
Case of Small Delay Residual 



-q 

co 


RUN 159-21-42 NU. RcC.= 5 i.HtC.NO.e= 1 

50UKC5 ATE-3 ?-R NO. PT3.= *•'' S. RtC. NO. »= 1 

SLUM. JELAY = 411 4. 2' 1 - 

INST, CL LAY = -5531: 


OAV= 1E9 T iMA = £i.1> ••• 22. 1.1412 


F Kt£0 

1 9.90"?'' "P 

* {?'*■•.> r. 

3 

4 l^!.2^">rr o« ; 

5 K*.3'C v-nr. 

6 i n. 4 rrf 

7 i f .s^rrv. 

a 1 

9 1 0.7"'Cf‘'' A 
i r i f •a:, r co^i} 

11 1C .9 n 0C v 

12 1 1 ,rrf rf.r»n 

13 11 .I'^r^r.*) 

14 1 1 

15 II 

16 U.4^.rVl 

17 1 1 ,5™''-T' r 

18 n.6rcrfl.^ 

19 11.7!'CC0*C 
2* ll.a'^Ot*© 

21 11.90"'* CO 

22 12 CO'*" 

23 1 2.1 or oo '-o 

24 12.2" rccc 
2b 12 

26 1 2 • 4 * - * -* •: n 

27 12. S*"' 1 '': 

2 8 1 2 .6 '".C-' 'T 

29 12.7 

30 12.8f"' , t'f' * 

31 12 .9 K-'" 

3 2 13.0*?f *“">1 

33 1 3. IOC .*>0"C 

34 1 2 m 2 

35 12. "COO" 

36 12.4 ''CCnf*.:' 

37 

38 13.6 ' ^ *V*. 

39 13.7 

4 *5 1 3.0v .CO';" 


AMPLITwt 
' • C 3^ 54 7461 S' 4 
•' • *>55 3* 11 >27350 
'■> • 13 I4fc2,.25^6 ! > 
'•• 1 $9 12 It* 67c 39 
•*•2 5 "6 494 96226 
28; 71 7701*313 
C .2e3'3*67€235 
r ’ • 252 2* 525 r < 92 
lS'S*' 652^36 
n. 11 '647^45-ji7 
i rr £f 166031 7 
<'• . 231 o£. 8966706 
r . 29 94 £6736243 
•^.3>J '‘24o4 59 Oil 
A . 761 37 1 ?>fc 3 1 1 6 

S3 l*<ni 69269 

1 . '*S13"6'* S44 1 

1 • 2^5'>442i 97"' 9 

1.294355521 *96 
1.241 732 34 299 , ’» 
1.343 tco9993'"l 
1 . 2^2=61212^26 
1 . 223379329313 

1 . 1 ' u " C 4 65 d 666 

9 5 *.53 971 26 '4 
0.76" 852 '21 232 
^.55 7 ! 49^52 155 
A • 4" 9336* 64646 
0. 23^3 •■'•5354 737 
^.<■75427^25122 
-. I'' 5 3 3C^ r 3589 
'.226 7 ■' 36 3 o 2 6 4 
6 . 31 - 7554676 1 4 

3c .! l£< 7C52’'l 
. 3 7 3 2 4 5 3 54 64 9 

■>.31 J6 667 <*752 5 
« .267514605038 
”■ .163625 5 26 151 
' • 1 *45(j'S3ei 


TAU 

- 1 • 55561 9 1 o24 ! r ' 
-3. 33371 7969549 

- 1 • 4l 371 3 1 oof-69 

- 1 . *♦62 1 36 >2 4 o76 
-1 • 4751 573'* 1 691 

- i * 461 753 70 64 7i 
-l . 46o6969o r 197 
-1. *33473 1 73477 
-1 . 45o834 1 75527 

- 1 • 56 '*48994 3651 

-1.461 6«*i 1 1 o9o7 
-1. 5^ ">3 35 o ‘•4552 
“1. 3*4322 6'-7V72 
- 1 • 3 "• 1 5 **3 7t 7 o 1 5 

- 1 • 5 n 5 £.21 1 76932 
-1.5159565 *>6433 
-1. 51 *.l‘'997boo3 
-1 • 497046 735695 
-1 • 49 524 '• 34* *2 69 

- 1 • 53 7604 1966 70 

- 1. 51476312067*1 

- 1 . 5*7552977^60 

“ 1 . 45733 5**co7 97 

- 1 • 4 5 365** 33 o 79 2 

- 1 . 492 J'>bl 7 67 30 

- 1. 48*..- 76<*1 7747 

- 1 . 4 5 j 0 99 /r, 9 3 066 

- '. 54-,7l t 53479 

- 1 .55717** 52 "ie. 64 
- i • 04812 1 7 o97 79 
“1.53 "o? b 1 6 -'646 
-1 • b 2 4 r - ■> 3346123 

1 • 61 61 42 r 1 792 3 
1. 512675 *4 5^:01 
1 • 5 r 7 5 4 9 f d S o 3 7 
1 . 4995 o663joo 6 
1 • 4ciofc4cl 57o9 
1 • 433069536573 


C0WK2L AmPL 
r > .^1 2l 95 416460 I* 

C . '152 5oc4 > ’»3o 6 Iv 
0. ft 41 d'-'66‘*4 378 I 
! .'0459991 16163 I 
r . '*3189 142441 1 I 
v. 051 523 31 7756 1 
' <‘•*91 672655 2:* 7 I 
2- • '61826 324793 l 
A .«- 61 1 2 1 606763 I 
'3174b 16912* I * 
«"• 023926691497 I v 

* • '662959 r 9303 1 

' . 122846953331 I 
*. 15*9 3 J 759923 1 
'• • 2 3857*7360<3 2 1 

*■ • 2'-. 1 3 6- 305276 I 

0.337241723422 I 
C .87 04 33 9 ■' •"* 1 3 . • 1 

* • 4" f. '49 350271 1 

' • 4 2 ■'*• 2 1 9 7 1 7 1 **7 1 
"•Mat 41 17*34*' I 
*’•4'' 51 '*-3 7628c 6 1 

•'♦ 365534 "»13574 I 
* • 35 ' 7 1 1 v S47 1 6 I 
*' .33 57 8922961 6 I 
‘ . 25.J^3**3624l o 1 
' . I'?cj58(j6£32b r> I 
'• 1 394754 702^.4 I 
".^844797996:7 1 

0. ' , 3odo7o 76'19 l 
‘•.■'3625125"' '-63 1 
f% . " 6 1 4 7b3o97 2 1 1 

r .0791 3’' 872: J8 I 
> . , '94'' 573 n 48ii 1 

*'.'•9 97 44 2 37313 1 
'• '9655223821 * I 
’ . *'• d ‘5 9 1 2 9 c i 3 1 9 1 
C .t»o9 5 l'*229o42 1 


'• . ' 49o5778324 9 I 
••.02627515670 3 I * 



Figure 20c. Fringe Spectrum for Pseudo Random Sequence Modulation 
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Figure 21a. Correlation Function for Sidetone Ensemble Modulation - Case of Large Delay Residual 
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Figure 21b. Power Spectrum for Sidetone Ensemble Modulation - Case of Large Delay Residual 
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Figure 22a. Correlation Function for Sidetone Ensemble Modulation - Case of Small De lay Residual 
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Figure 22c. Fringe Spectrum for Sidetone Ensemble Modulation 
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Figure 22c. Continued 



central correlation peak with small side peaks is a characteristic of this 
modulation. Figure 19a is the correlation function of the 1 Mbit/sec rate PRS. 
This modulation exhibits a very broad peak with negligible side peaks. The 
expected bias towards higher frequencies can be seen in both the spectra (Fig- 
ures 18b and 19b). Figures 20a, b and c correspond to the case of PRS modula- 
tion where the bit offset is once again less than one. However, the typical fea- 
tures of the correlation function and power spectrum are identical to those of 
Figures 18 and 19. Figures 21 and 22 are the results of a sidetone ensemble 
run. The range tones are clearly visible as can be seen in the spectral plots. 
The correlation function exhibits the strong side lobes that are characteristic 
of this modulation. 

Analysis of the characteristics of the five types of modulation show that 
broadband noise, either the spacecraft front end noise (this is the noise generate 
by the transponder front end when it is turned on; the transponder is not 
triggered by any ground emitted signal) or the transponded noise spectrum whicl 
is originated on ground is the most advantageous type of modulation. This is 
predictable theoretically since a broadband noise signal contains the most 
frequencies within a given passband, and one expects this to make the central 
fringe stand out most prominently. In practice, this was confirmed, first by 
the fact that the fringe phase slopes for these two types consistently showed the 
least amount of jitter and secondly that the correlation peaks in these two types 
stand out strongly from the background. 

t In our data processing scheme, a standard integration schedule of 5 
records was followed. Each digital tape has approximately 3 minutes or 900 
records. In cross-correlating a pair of tapes, averaging is done over 5 records 
and the whole tape is processed in sequence of 5 records. Thus one has a delay 
and delay rate parameter approximately every second (5 records). The delay 
and delay-rate (fringe-rate or differential doppler) are printed out and plotted 
with all the accessory information. Figures 23,24, 25, and 26 show samples 
of the processed runs. The data (time-delay and fringe-rate) shown in these 
runs are smoothed and used as inputs to the orbit determination program. 

Table 9 shows a sample list of satellite runs that were processed. The 
time delay and the doppler correspond to the time of the run. The format of 
the run indicates the day of the year, hours, minutes and seconds of universal 
time. Also, indicated in the table are the various signal modulations corres- 
ponding to the data run. 


79 



00 

o 


RUN 150-21-' 7 

SOUPCF ATS-t Fi-X 


POSMAN - M 3 J A VE 
NRFC = 5 


DAY 

153 
t 58 
1 58 
1 58 
1 58 
1 58 
1 56 
158 
1 58 
1 58 
158 
156 
156 
158 
158 
158 
1 58 
158 
158 
l 58 
1 58 
158 
1 58 
158 
158 
158 
1 58 
1 58 
1 58 
156 
158 
158 
158 
1 58 
1 58 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
153 
158 
1 58 
I 5R 
158 


HP MIN 
2 1 7 


2 1 
2 l 
2 1 
21 
2 1 
2 1 
21 
2 1 
21 
21 
2 1 
2 1 
2 1 
21 
2 1 
21 
2 1 
21 
21 
21 
21 
2 1 
2 l 
21 
2 1 
21 
21 
21 
21 
21 
21 
21 
21 
21 
21 
21 
21 
2 1 
2 l 
21 
21 
21 
2 1 
21 
2 l 
2 1 
21 
21 
21 


StCONP 

1 .3 86 .’ < 

2. 39 4v 
3 . 4 -.'2 

A. 4 l-:.'. o 

5.4 1 8< f 

6.426 ( 

7.4 34 0 

8 .4429 r 

9 . 45 c r ( 

l«; .458. < 

1 1 .4 66< 0 

12.4 74—1 

13.4 Q2«‘ * 
1 4.490 r 
1 5. 4 98' «' 
16.5C6C 0 
17.51399 
1 8.52199 
1 9.53'JOC 
2'’' .53799 
21 .54599 
22.55399 
23.562- 0 
24.56999 
25.57799 
26.586 -»£ 
27.59399 
28.67199 
29.69999 
3 .6 18'* 
31 .62599 
32.63399 
33.642*0 
34.64999 
35.65799 
36.66599 
37.674C 0 
38. 681 99 
39.68999 
46 . 6 98' C 
41 .70599 
42.71399 
43.72198 
44. 7 3 0^ 0 
45. 73799 
46.74599 
47. 754C 9 
48.76199 
49.76999 
50. 7789 C 


DOPPLFR 
-23.275 
-23.210 
-23.225 
-23. 42C 
-2 3. 1 8C 
-23.230 
-23.390 
-23.155 
-23. 260 
-23. 265 
-23.225 
-23.24C 
-23.230 
-23.225 
-23.165 
-23. 25C 
-23. 260 
-23.25* 
-23.175 
-23. 1 75 
-23. 210 
-23.060 
-23. 355 
-23. 160 
-2 2. 205 
-23. 220 
-23* 1 80 
-23.215 
-23. 230 
-23.150 
-23. 235 
-23.2CC 
-23. 245 
-23. 235 
-23.21* 
-23. 165 
-23.095 
-23. 28C 
-23.1 70 
-23.1 9C 
-23. 295 
-23. 285 
- 2 3 • 1 4C 
-23. 245 
-23. 225 
-23.1 e* 
-23.245 
-23.200 
-23.215 
-23.250 


■3250 
■82 50 


T AU 

825<* .467 14 ^ I 
82 5 f ' .4 3262 5! 
62 50. 5*4016 I 

44 8'- 57 I 
5720231 
- 82 50 . 49 2 1411 
-8250.3461 121 
-8250.413295 I 

— 82 50.36557 5 1 
-8250.4310551 
-8250.4549221 
-825u . 430 142! 
-8250.3792651 
-825c .3688721 
-8250.423093 ! 
-82 50.382 74 4 I 
-8250.4163661 
-825C .3952031 
-8250 • 3896461 

— 825*' .3341621 
-8250.3401991 
-82 5C .366431 I 
-8250.349954 I 
-825C. 349 74 1 I 
-8250.3350991 
-825C . 3288671 
-8250.291723 I 
-8250. 3119331 
-325C.314 10 31 
-3250.3405451 
-8250.3138731 
-8250 .2870361 

— 82 50 • 30 3099 I 

— 82 5T.274545I 

— 825 0.2357891 
-8250.2659431 
-8250.2744531 
-8250 .263706 I 
-825C .2508161 
-8250.21 12251 
-82 59 . 258621 I 
-8250.2621431 
-825C .279181 I 
-8250.2248391 
-8250 . 1814461 
-8250. 2203121 
-8250.2205331 
-8250.227134 I 
-8250.2224751 
-8250 .1961631 



VLB 1105 VAR 1C 5/VAM 105 SAT-F 
TAU 


* 

* 


* 

* 


* 


* 

* 

* 


* 

* 


4 


* 



Figure 23. Fringe-Rate and Time-Delay Versus Time for the ATS-3 Satellite - Noise Modulation - 
Run No. 158-21-07. Rosman - Mojave Experiment . (Sheet 1 of 4) 
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Figure 24. Fringe-Rate and Time-Delay Versus Time for the ATS-3 Satellite - Spacecraft Noise - 
Run No. 158-21-54. Rosman-Mojave Experiment. (Sheet 1 of 4) 
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TABLE 9 

ATS-3 VLBI Runs Processed* 


Run^ 

Signal- 

Modulation 

Delay 
( Sec) 

Doppler 

(HZ) 

152-19-38-31 

Noise Translation 

4076.454 

18.460 

152-19-39-11 

Noise Translation 

4076.242 

18.270 

152-19-39-58 

Noise Translation 

4076.075 

18.310 

152-21-49-00 

Noise Translation 

4102.998 

11.910 

152-21-50-30 

Noise Translation 

4102.745 

11.390 

152-21-51-59 

Noise Translation 

4102.495 

11 .450 

153-21-37-01 

Spacecraft Noise 

4118.150 

12.235 

153-21-38-30 

Spacecraft Noise 

4117.898 

12.210 

153-21-40-00 

Spacecraft Noise 

4117.623 

12.010 

153-22-39-31 

Noise Translation 

4113.592 

7.655 

153-22-41-00 

Noise Translation 

4113.420 

7.565 

153-22-42-30 

Noise Translation 

4113.259 

7.425 

153-22-52-01 

Pseudo- Random Sequence 
(500 KB) 

4110.170 

6.665 

153-22-53-30 

Pseudo-Random Sequence 
(500 KB) 

4110.144 

6.215 

153-22-55-00 

Pseudo-Random Sequence 
(500 KB) 

4110.058 

6.460 

154.22-33-00 

Spacecraft Noise 

4112.925 

7.925 

154-22-33-45 

Spacecraft Noise 

4112.847 

7.950 

154-22-34-28 

Spacecraft Noise 

4112.812 

7.735 

154-23-37-01 

Noise Translation 

4114.256 

2.440 

154-23-38-20 

Noise Translation 

4114.197 

2.330 

154-23-39-40 

Noise Translation 

4114.140 

2.130 

158-21-55-01 

Spacecraft Noise 

4115.624 

9.600 

158-22-06-31 

Noise Translation 

4116.889 

8.710 


*Each run is of 9 minutes duration 
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TABLE 9 (cont. ) 


Run * 

Signal - Modulation 

Delay 

(sec) 

Doppler 

(Hz) 

158-22-29-31 

Pseudo- Random Sequence (1MB) 

4119.491 

6.940 

158-23-15-49 

Noise Translation 

4121.046 

2.990 

159-12-40-31 

Spacecraft Noise 

3985.508 

12.710 

159-13-35-01 

Noise Translation 

3999.886 

12.345 

159-15-29-31 

Noise Translation 

4020.230 

18.990 

159-15-41-31 

Spacecraft Noise 

4023.417 

19.380 

159-16-33-00 

Noise Translation 

4042.300 

20.600 

159-18-37-00 

Spacecraft Noise 

4076.495 

19.760 

159-19-32-31 

Noise Translation 

4092.300 

17.590 

159-21-11-22 

Noise Translation 

4110.034 

12.630 

159-21-22-01 

Pseudo- Random Sequence 
(500 KB) 

4112.733 

11.830 

159-21-33-31 

Pseudo-Random Sequence (1MB) 

4114.183 

11.045 

159-21-45-01 

Spacecraft Noise 

4115.940 

10.170 

159-21-56-31 

Sidetone Ensemble 

4118.090 

9.195 

159-22-11-31 

Noise Translation 

4119.702 

7.720 

159-22-27-01 

Spacecraft Noise 

4121 .421 

6.760 

159-22-39-31 

Pseudo-Random Sequence 
(500 KB) 

4122.419 

5.650 
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6 . 


DETERMINATION OF THE ATS-3 SATELLITE ORBIT PARAMETERS 


Two approaches have been followed to determine the orbital elements 
of the ATS-3 spacecraft using the VLBI data. The first one is an analysis based 
on the first order orbit theory that is developed exclusively to fit the VLBI data. 
The second approach is the utilization of a highly sophisticated orbit determina- 
tion program - the "GEODYN" program. Both the approaches and the corre- 
sponding results are presented in detail. 

6.1 ORBIT ANALYSIS 


6.1.1 Formulation of the Geometrical Delay 


The geometrical relationships involved in the problem are shown in 
Figure 27 . The satellite is located at S relative to the Earth's center, and at 
R, and R 2 from^ground^ stations A and B, respectively. Stations A and B in turn 
are located at P 1 and P 2 relative to the Earth's center. The baseline D is simply 
D = p i “ p 2 * The geometrical delay, At, is given by the equation 


D R 2 Ri • R-2 

? At = + RT + -fC—- R i 


DS DP., 


Ri 


Ri 


Rj • R 2 


- R, 


(43) 


For | S |-»°°, i.e. an infinitely distant source, Equation 43 reduces to cAr = D-S. 
This is the familiar form of the delay equation used injstellar VLBI. For |S| 
large in comparison with the Earth's radius the term P_Jf i s still strongly dom- 
inant in Equation 43, and the remaining terms can be viewed as "corrections." 

D • P * 

the factor -j— ^ representing the contribution of an eccentric baseline vector 

and the expression within brackets accounting for the finite height of the satellite 
source. 


If Lj denotes the longitude of station i from the satellite's ascending 
node (£2) and X. its latitude, then 


D = (cos Xj cos Lj - cos X 2 cos L 2 , cos Xj sin Lj - cos X 2 sin L 2 , sin Xj - sin X 2 ) 


(44) 


In the following analysis, we assume the Earth's radius to be unity, and all 
distances are normalized to it. 
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Let to e represent the Earth's angular rotation rate, anda ; the longitude 
of station i measured from the ascending node at the time of perigee passage, 
which is taken to be the epoch in the present analysis. Thus L i =GJ c t + a , ; * = 2 - 

Substituting these expressions in Equation 44, and expanding and collecting terms, 
one obtains 


D = (A cos [co e t + a] , A sin [co e t + a] , sin Xj 


sin X 2 ) 


(45) 


where 


A = [(cos X t - cos X 2 ) 2 + 4 sin 2 



cos x 2 ]^ 

(46) 


and 


tan a = 


cos Xj sin Oj - cos X 2 sin a 2 
cos Xj cos aj - cos X 2 cos a 2 


(47) 


From Figure 28 it can be seen that (sin X t - sin X 2 ) is the length of the projection 
of D onto the polar axis; A represents the length of its equatorial projection; and 
a is the direction of this equatorial projection relative to the line of nodes at 
t = 0. One can describe S, the satellite's position, in terms of u, the argument 
of the latitude and i orbital inclination as follows: 


S = (S cos u, S sin u cos i, S sin u sin i) 


(48) 


From Equations 45 and 48, the first term on the right side of Equation 43 
becomes 

^ [cos u cos (w e t + a) + sin u sin (co e t + a) cos i] 

R 2 R 2 s 

+ sin u sin i (sin Xj - sin X 2 ) (49) 


This is the explicit form of the first term on the right side of Equation 43. 
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Figure 28. Projections of the VLBI Baseline 

The first term on the right side of Equation 49 represents the contribu- 
tion of the equatorial component of the baseline D; when a y = a 2 , A = 0 and this 
term vanishes. Similarly, the second term on the right represents the contribu- 
tion of the polar component of D; when \ = \ 2 , it vanishes. 

The second term on the right side of Equation 43 is - , and its 

numerator is a constant that depends only on the locations of the two stations. 
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The third term can be written as (Reference Figure 27) 



D • R 

where cos 6 = - — - . This term will be treated later. 

DRg 

6.1.2 The Synchronous Satellite Orbit 


For the synchronous satellite orbit we assume i 2 « 1 and e^c< 1, from 
which it follows that sin i ~ i and cos i ~ 1. For small eccentricity, we can 
approximate u, the argument of the latitude, in Equation 48 by the well-known 
form: 

u = cj = £2' t + 2e sin S2' t 

(51) 


where oj is the argument of the perigee (measured from the node),fi' is the 
satellite's mean motion, and t is the time measured from perigee passage. 
Expression 49 becomes 

^ ^ [cos (oj t + a)] + — isinu (sin Xj - sin X 2 ) 

R 2 R 2 R 2 

Applying Equation 51 and expanding the cosine, we get 


D • S AS 

: -jr— [cos (a - a>) + rjt sin (a - co) + 2e sin (a - co) sin co e t] 


R 2 R 2 


s 

+ i sin u (sin X, - sin X,) 

r 2 


(52) 


This is for the case of a nearly synchronous satellite orbit. The symbol^ 
represents ' - co e ■ 

6.1.3 The Rosman - Mojave Baseline 


The coordinates of the stations are: 

Longitude Latitude 

277J124 E 
243?112 E ; 


Rosman: 

Mojave: 
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35H99 N 
35°332 N 



Latitude difference is 0.132°, i.e. the baseline lies in a plane essentially 
parallel to the equator. We shall, accordingly, ignore the term as well 

as (sin Xj - sin X 2 ) and (cos X 1 - cos X 2 ). Then the last term in Equation 52 van- 
ishes, as does the first term in brackets within Equation 46. This latter expres- 
sion becomes 


d 1 ^ 

A = 2 sin(— -)cos Xj = D, 


while Equation 47 reduces to 


sin a. - sin a , 

tan a = 

cos aj - cos a 2 



Equation 52 reduces to 



[cos P + rjt sin P + 2e sin P sin u> e t] 


(53) 


where P denotes the difference ( a- co). However, the nodal regression and 
apsidal rotation of the satellite's orbit are both first order effects in J 2 , and P 
should, therefore, be viewed as P = P Q + P 0 t, where P Q represents the sum of 
nodal and apsidal rates. Equation 53 becomes: 


PS 

Rt 


A [cos P Q + (r? - P Q )t sin P Q + 2e sin P Q sin w e t] (54) 


The first term on the right side of Equation 43, therefore, takes a 
simple form. The third term can also be simplified accordingly. Equation 50 
can be expanded in powers of the small parameter — < i , and the series 
truncated after the first two terms: 1 


R 1 R 2 D 2 

-nr _R i * 


(55) 


Since 


cos 6 = 


D ■ R 2 


D(S-P 2 ) 

dr 2 


D S 
DR-, 


+ (higher order terms) 
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one can write 


R i R 2 

R, 


D z 


- R 1 = 


2R, 


■ 


\DR, 


cos 2 P_ 


+ higher order terms 


(56) 


D z 

2Rj sin 2 P Q 


to first order in small quantities . 

The following expression for the geometrical delay is obtained when 
equations 54 and 56 are inserted into equation 43: 


-T = 


o/ r 2 \ 
cos p o -Ts\i rj 


D P 0 

sin2 p o - "5s - 


(57) 


+ (tj - P Q )t sin P Q + 2e sin P Q sin co e t 


where T denotes the expression 

6.1.4 Extracting the Orbital Elements 


cR-, c R 0 

— - — — . — . At 
AS D S 


A best-fit diurnal curve is obtained to fit the observational values of 
time-delay and fringe-rate using the least squares technique. In accordance 
with Equation 57, the form of this curve is chosen to be 

T = T o + Tj co g t + T 2 sin (o^t' + T 3 ) (58) 

where T Q> T Jt T 2 , and T 3 axe determined by the curve-fitting procedure. The 
variable t' represents time as measured from some convenient reference, such 
as 0 h UT on a specific day. When co e t’ = -T 3 , the satellite is passing through 
perigee, and this value of t' (viz. t’ = - T 3 /uj is the epoch adopted above. 

Equating the expressions 57 and 58 and evaluating them at t = 0, as 
well as their derivatives at t = n / 2 and t = 0, one gets 


T = cos P o - K sin P - M 

0 o o 

(v-p o ) 

T, = — ^ P 

1 ° 


T 0 = 2e sin P 

L O 


(59) 

(60) 
( 61 ) 
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where 


and M = 


K = 


D R 2 
2S RT 


d-p 2 

DS 


Equation 59 is a quadratic form in cos P whose solution is 


P o = hk L 


-1 + / 1 + 4K(T q + K + M) 


( 62 ) 


WithP known, both 17 - P 0 and e can be found from Equations 60 and 61, viz: 


v-P n = 


w e T l 


0 sin P„ 


(63) 


e = 


2 sin 


(64) 


and, as discussed above, 


t = 


CO 


(65) 


Equations 62-65 constitute the first-order solution for the elements 
of the satellite's motion. 

6.1.5 Determination of ATS-3 Orbit Parameters 


From each tape pair of VLBI data, a representative set of time-delay 
and fringe-rate parameters are tabulated. A curve of the form of Equation 58 was 
fit to the time-delay values obtained, using least squares. The curve is shown 
in Figure 29. If we denote by A r i , the quantity 2. l.j. for i = 0 , 1, 2, and 
At 3 = T 3 , then from the curve as shown: c 2 

Ar Q = 4039.7 ^sec 

Atj = .005 gsec/deg 

A t 2 = 71.32 jxsec 

Ar 3 = 91.308 deg 


104 



DELAY (mj) 


ATS -3 DAYS 152-154. 158-159. 

ZERO TIME= ZERO HR. DAY 152 
DELAY= A + B» FR • T+ C • SIN (FR • T + D) 
FR- 361.0 DEG/ DA 
A = 4039.6840 MiCROSEC 



Figure 29. Time-Delay Curve for ATS-3. Days 152-159, 1971 


The calculated value for ^ is 

T q = .3595 

T i = .487X10' 6 (deg )’ 1 

T = .0063 

2 

T 3 = 91.308 deg 

D • p. 

For the Rosman-Mojave baseline, D s = - 0.036. Inserting the above values 
into Equations 62 through 65 we find: 

P Q = 66°20' 

t?-P q = .011 

e = .0034 

t' = 18 h 00 m 

The term 77 = fi' - <u e has a periodic variation that reflects the variation 
in £ 2 ’ (the satellite's mean motion) induced by dynamical resonances at the 
synchronous altitude. The official ATS-3 Project orbital parameters, for 
instance, cite values of 1435.98 and 1436.125 min. for the anomalistic period on 
8 May 1971 and 31 May 1971, respectively. Thus, the least squares fit will 
"average out" this variation in rj and the above equation for 17 -P 0 must be 
interpreted as an equation for -P Q , i.e. the sum of nodal plus apsidal rates only, 
rather than for tj-£ 0 . In the limit, as i-»0, = -2 £2 = +2 a, and that (tb-a) = 


Thus 
S 2 = -.001 

If we restrict the delay values to those obtained during a 24-hour period only, 
then rj will not be averaged by the least squares fit, but rather will assume its 
"instantaneous" value for that day. Accordingly, we fit a curve to the data for 
days 158 and 159 only. This yields a value of At = .3767X10’ 3 , which is appro- 
priate to about noon on day 159. Then 

Tj -P 0 .3767X10' 3 

— = = .0685 

-P Q 5.5X10 - 3 
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from which 77 = -.0102 deg/day. If T and T' denote the sidereal day and the 
instantaneous anomalistic period of the satellite (both in minutes of UT), then 
T' - T = -15 17 = .153. Thus T' = T + .153 min., where T = 1436™068. Finally: 

T' = 1436 T22 

In conclusion, table 10 compares values for the ATS-3 orbit elements supplied 
by the Project Office, against the first-order values obtained herein: 

TABLE 10 

Comparison of ATS-3 Orbital Elements 


Orbit Element 

VLBI Value 

ATS Project Value 

Period (day 151) 
Semi-major axis 
Eccentricity 
Time of perigee 
Position of perigee 
Nodal regression rate 

1436 m .22 
42168 km 
0.0034 
1 8 h 00 m 
0.5 deg 

-0.011 deg/day 

1436 m . 125 
42165 km 
0.0029 
1 8 h 05 m 
359.288 deg 
-0.0134 deg/day 


6.1.6 Error Analysis 


Errors in the measured values of relative delay and delay rate arise 
chiefly from two sources, tropospheric perturbations and equipment misalign- 
ments., These effects manifest themselves as errors in the A t { (i = 0, 1, 2, 3) 
that are determined from the least-squares fit of Equation 58 to the delay rate. 

Equations 62 - 65 can be used to express errors in the orbital elements 
in terms of errors in the Ar ; . For simplicity of exposition, we shall temporarily 
replace Equation 62 by 

T 0 = cosP o ( 62 ) 

Then, as is readily apparent by taking derivatives of Equations 62' and 63-65 
with respect to the : 
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CL? 

0 

1 

cot P Q 

3Ar c 

1 

> 

O 


3P 

P cot 2 P 

3P 

P 

O _ 

O 0 

0 

0 

3Ar 

0 

At 

0 

’ dA Tl ~ 

Ar l 

3e 

e cot 2 P 

0 

be 

e 

3At 

0 

At o 

’ 3Ar 2 ~ 

At 2 


3t' 

t' 

3A t 3 

cn 

I-* 

< 


Clock-epoch misalignments produce a constant offset in A r 0 , and their effects 
on P 0 , P Q and e are analytically described by the expressions for d(P 0 ,p 0 ,e) 
cited above. Differential atmospheric path lengths may also produce constant- 
bias errors in Ar 0 and these will contribute additional errors to Ar 0 . Small 
local oscillator (L.O.) offsets between pairs of terminals will lead to an error 
in Ar , and this will have an effect on the derived value of P Q , as described 
above. 

Consider now any errors in the adopted baseline D. These arise 
because the positions of the antennas are not precisely known. The constraint 
that X 2 - X 1 limits the consideration of positioned errors to those lying in the 
plane parallel to the equator at latitude . Thus, only the difference in longi- 
tude between the two stations enters into the analysis. The baseline error, 
being one-dimensional, can be equivalently expressed as an error in its length, 
D. From the definition of T. , we have ELEi T| . Again employing Equations 

62, 63 and 64: * 3d D 


3P„ cot P„ 
o_ _ o_ 

3D D 


3D 


D sin 2 P„ 


3e _ e 
90 Dsin 2 P 0 
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6.1.7 Second Order Effects in the Orbit- Analysis 

The orbit analysis is extended to include second-order terms in the 
small quantities and ~ , where D is the baseline length, and a is the satellite's 
semi-major axis (= 6.61 earth radii), and P 2 is the distance from the Earth's 
center to station 2 (= one earth radius). As before, e 2 and i 2 are neglected. 



Figure 30. Geometrical Relationships - Second Order Effects 

For the general non- equatorial baseline, equations 43, 50 and 52 may 
be combined to yield: 

s • 

CAT = T^i A t C0S P C + 1 Sin P o + 2e Sin P o Sin W e t] 

dp 2 

+ i(sin Xj - sin X 2 ) sin u + — — (1 + e cos [u w])| (66) 

D , 

+ 2R; sin d 
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2 


If \p the angle between P 2 and S, 


from which 


(f) ■ 1 - 2 (t) cos * + ("f) 

- 1 + (f) ^ * 4(t) -4(f) cos ** 


(67) 


where (Ref. Figure 30), 

cos \p 


= cos v cos u + sin u cos f 
= cos X cos(L 2 - u) + i sin X sin u 


( 68 ) 


and 


p 2 p 2 • 

g-= — 0 + e cos [u - co]) 


(69) 


Inserting equations 68 and 69 into equation 67, and reducing the terms that 
result: 

— = 1 + |cos X cos(L 2 -u) + e cos (u - co) cos (L 2 - u) cos X + i sin X sin u j 


cos X cos (L 2 - u) 


Equation 70 can be inserted into Equation 66. 

Returning to the last term in Equation 66, and using Equation 55: 

/D S \ 2 /D S \ /D P 2 \ 

COs2 6 = \DR 2 / " 2 V DR^/ \ DR 2 / 

Employing equation 54 and simplifying: 

= 1 j | cos 2 P Q + 2 jcos 2 P Q cos X cos(L 2 - u) + (t? - P Q ) t sin 2P Q 


sin 2 d 


+ 2e sin 2 sin cjt 
o e 




cos P„ 
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Finally, the term can be written as 


D D / s \ D f 
2Rj = 2S \R^ / 2a [ 


1 + e cos (u - go) + 2 — cos v 
a 


( 71 ) 


where v is the angle between and S, in Figure 27 . 


Equations 67-71 can be inserted into Equation 66 to obtain an explicit expression 
for cAr which is accurate to second order in the quantities cited above. Equa- 
tion 66 supersedes Equation 57. To first order, it is: 

( DP, 


cA t = < A cos P Q - 


2 » 1 -(£) cos2p <>] + A (r) 


COS P Q COS ?2 COS X > 


A(t? -P Q )t sin P Q 


2Ae sin P Q + i(sin \ - sin X 2 ) cos w 
2Ae cos P Q + i(sin Xj - sin X 2 ) sin co 


sin Wgt 


cos co e t 


(72) 


The second order expression contains constant terms, as well as first and 
second harmonics in co e t. The constant and linear terms are ignored here 
because they add little of qualitative interest to the solution. The terms 
involving sin co e t and cos co e t contribute an important correction to e. They are 


+ ^] cosw e t + e [ 2 A (” ) cos X sin (P Q + P 2 ) + a ( D ) sin 2P 0 | 
sin co e t + i(sinXj -sinX 2 ) jcos co sin w £ t + sin co cos co e tj 


(73) 


For the Mojave-Rosman baseline, the last tertn of expression 73 is negligible 
and the remainder is of order 


-.1 De cos oo e t + .2 De sin co e t 
To second order, then, Equation 61 should read: 

g At = 1 2.2 sin a) g t + 0.1 cos a> e tj e = 2.2 sin (co e t + 2.5°) 
from which the corrected value for e would be 


e = .0031 

which compares more favorably with the value of .0029 given by the ATS 
Project than the value of .0034 as determined by first-order theory. 
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Also, the value for T 3 obtained by the least-squares fit is 91.2 deg. At the time 
of perigee then, 2.5° = fl t'^ +91.2°, where t' is time measured from O UT. Thus, 
the value of UT at perigee passage is t^ = -5 n 55 m . In other words, the time of 
perigee passage is 18 h 05 m which agrees exactly with the figure supplied by the 
ATS Project Office. 

The second-order terms involving the second harmonics are 

A /P 2 \ 

ei (sin Xj - sin \ 2 ) sin (2co g t + co) + — I — I e cos P Q cos (2co e t + co-L 2 ) (74) 

For the Rosman-Mojave baseline the first term in expression 74 is negligible 
in comparison with the second. The second term is of order 

10 4 D sin (2co e t + 35°) 

By extending the least-squares fit to the data for days 152-158 to include a 
second-harmonic term of the form b sin (2co e t + B), this term takes the form 

2 x 10‘ 4 sin (2oo e t + 4°) 

which indicates reasonably good agreement in the magnitude, but rather poor 
comparison in phase due to lack of sufficient observations at different phases 
of the satellite orbit. 

6.2 ORBIT DETERMINATION USING THE "GEODYN" PROGRAM 

This program is a sophisticated orbit determination program which 
has been modified to accept the satellite VLBI parameters (time-delay and 
fringe- rate) from an arbitrary number of baselines as inputs to the program. 

The GEODYN program has the capability of estimating that set of 

• Orbital elements/determination of definitive orbits 

• Station positions/ geodetic parameter estimation 

• Measured biases/tracking instrument calibration 

• Force model parameters 

which give the best fit for orbital tracking data. The program is operational on 
the GSFC IBM 360-95 computer. The total number of adjustable parameters 
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affecting any one arc may be as large as 70. The total number of adjustable 
force parameters affecting any one arc may be as large as 20. 

The Inputs to the program are: 

• VLBI Parameters (time-delay, r g and fringe-rate v f ) 

• Measurement biases 

• Station coordinates 

• a priori satellite cartesian elements and epoch time 

• Geopotential coefficients through degree 22 and order 22 

• Solar-flux and geomagnetic data 

• Satellite drag-coefficient, satellite emissivity and mass 
The outputs from the program are: 

• Double precision adjusted cartesian elements; correlation 
coefficients between elements; standard deviation of fit 

• Kepler elements 

• Adjusted geopotential coefficients 

• Adjusted Earth-fixed cartesian and geodetic station coordinates; 
correlation between determined coordinate values; standard 
deviation of fit for coordinate values 

• Instrumental biases - timing and constant biases 

• Summary of residual statistics 
Other features of the program are: 

• The motion of the polar axis is described by a mean pole position 
averaged over six-year periods since 1900; and by libration and 
periodic terms 

• Luni- solar perturbations on satellite motion 
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• Orbital perturbations due to earth tides 

• Precession and nutation effects of Luni-solar and planetary 
perturbations 

• Planetary perturbations on the plane of the ecliptic 

• Tropospheric and ionospheric refraction corrections 

• Solar radiation pressure 

The numerical integration technique that is used in the program is 
the "Cowell Predictor" technique. More details about this program can be 
found in the document titled "GEODYN SYSTEM, Vol I, II and III" by Wolf 
Research and Development Corporation, NASA, GSFC Contractor Report, 1971. 

The station coordinates used in the program are the following: 

Station Latitude Longitude Altitude 

Mojave ATS 35° 19* 53" .78 243° 06' 43" .95 887 (m) 

(MOJATS) 

RosmanATS 35° 11' 56" .05 277° 07’ 27" .30 828 (m) 

(ROSATS) 

The following table (Table 11) gives the orbital elements computed 
from the GEODYN program for various cases. It is understood that the orbital 
elements computed from the GEODYN program are more accurate. The orbital 
elements computed by using first order orbit-theory and a data-fitting procedures 
is crude and is presented here only to illustrate the VLBI orbit analysis. The 
last two columns show the orbital elements determined from two independent 
subsets of data and one can see a remarkable consistency in the results. 
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Epoch 71 y 149 d 00 h 00 m 00 s 

DATA- SOURCE 

ELEMENT 


TABLE 11 

■3 Orbital Elements from GEODYN Program 

Total Points 


843 for Time Delay 
843 for Fringe-Rate 


Semi-major axis (km) 

a 

Eccentricity 

e 

Inclination 

• ' 
i 

Position of Perigee 

(0 

Right of Ascension of 
the Ascending Node 

n 

Mean-Anomaly 

M 0 

Mean Orbital 
Longitude 

oj + Cl + M 0 

Bias Correction done 
by GEODYN 


Number of Data 
Points Used 


Number of Points 
Rejected 


Number of Iterations 
Period (mts) 


Nodal Regression 
Rate (deg/ day) 



Orbital 
Parameter 
from SHF 
R & RR Data 

Time-Delay 
{v g ) & Fringe- 
Rate (t'p) Data 
from All Days 

Data from 
Days 158 
& 159 
Only 

T n 9 v f 

Data from 

All Days with 

Large Residuals 

Rejected 

- Group 1 - 

42,165.43 
2.91 x 10" 3 

42,165.28 
2.879 x 10' 3 

42,165.17 
2.903 x 10“ 3 

42,165.270 
2.905 x 10" 3 

1 .706° 

1.818° 

1.712° 

1.712° 

359.288° 

2.620° 

0.223° 

0.277° 

82.214° 

78.563° 

81.223° 

81.195° 

87.150° 

87.443° 

87.172° 

87.159° 

168.652° 

168.628° 

168.618° 

168.632° 


No 

Yes 

Yes 


801 

578 

700 


— 

— 

72 


4 

5 

5 

1436.125 

1436.124 

1436.118 

1436.123 

-0.01340 

-0.01341 

-0.01342 

-0.01342 


2.905 x 10" 
1.712° 
0.277° 
81.196° 

87.159° 

168.633° 


1436.123 

-0.01342 
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7. DISCUSSION 


In order to determine the precision of time-delay measurements, 
each pair of data tapes was broken up into several subsets of data each ~ 1 
sec long (5 records) and the delay was computed for each subset separately. 
The distribution curve of time-delay from each pair of tapes was studied to 
determine the RMS Scatter. A total of one hundred and twenty pairs of tapes 
were processed in this manner. The delay-curve exhibited an average RMS 
Scatter of 5 nsec and the standard deviation of the mean is about 0.1 nsec. 

The precision of delay-rate measurement, dictated by the phase stability of 
the VLBI system, is estimated to be about 0.1 psec/sec in most cases. Since 
there was sufficient time-delay and fringe-rate data covering various portions 
of the orbit, the overall time-delay resolution can be estimated much more 
accurately. 

As described in an earlier section, the orbit of the ATS-3 spacecraft 
was computed using the GEODYN program. The GEODYN program insures a 
high degree of precision in orbit computation since it takes into account most 
of the conventional and non-conventional forces that influence the spacecraft 
motion. The orbital elements of the ATS-3 spacecraft were computed using 
various subsets of VLBI data. The semi-major axis is determined consistently 
to within 100 meters and the best value differs from that determined by using 
conventional SHF range and range-rate data by 160 meters. As shown in Table 
11, Group "A" and Group "B" (the last two columns) are the different subsets 
of data used in orbit computation using the GEODYN program. The agreement 
is remarkable. 


The arc- resolution (RAO ) may also be computed. By applying the 
usual interferometric equations, one can obtain 


RA d 


R cAt^ 
D sin d 


(75) 


where R is the range to the satellite, D is the baseline length, C is the velocity 
of light, Ar g is the time-delay resolution obtained in the interferometric 
measurements and q is the source-baseline angle. Using the following values 

R Q 

— — • — 12, A To =5x10 seconds and sin d = 0.5 

D 8 

RA0 — 35 meters 

For a time-delay resolution of 10 nanoseconds, the arc resolution is 
70 meters. 
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Similarly, the arc-rate resolution (RA 9 ) is obtained from the fringe 
rate measurements. 


= RA 8 


cAv ( R 
f sin 0 D 


(76) 


Where C is the velocity of light, Ap f is the fringe rate resolution ~ 100 p Hz 
and f is the signal frequency ~ 4 GHz 

RA0 ~ 0.1 mm/sec 

It can be summarized that the time-delay and the fringe rate measure- 
ment precision yield the ATS-3 projected spacecraft position within 100 meters. 
Higher precision can be accomplished if we use large recording bandwidths. 

We have the hardware/software needed to synthesize larger bandwidths and 
plan to utilize them in future experiments. With such a system, a time-delay 
resolution of the order of 0.1 nsec can be accomplished. An experiment in this 
direction has been proposed by us to the ATS C-2 project. 

The VLBI measurement of one way differential group-delay corresponds 
to the difference in the direct two way range delay measurements made 
simultaneously from Rosman and Mojave with SHF ranging transponders on 
board the satellite. In many ways, for equal accuracy of individual measure- 
ments, the direct-delays are more useful than the differential delays. However, 
the VLBI measurements yield very high angular precision. At present, VLBI 
differential group-delay data is about two orders of magnitude better than the 
conventional radio range data. When combined with conventional SHF range 
measurement, the VLBI parameters yield spacecraft position with very high 
precision. This leads to the discussion of the vector tracking concept in which 
three stations (two independent baselines) obtain two independent sets of time 
delay and fringe-rate by satellite observations. The two sets of VLBI para- 
meters are combined with a single range measurement from the master station 
(or the common station) to yield the spacecraft position vector with a resolution 
of the order of ten meters. 

VLBI measurement of fringe-rate or differential phase-delay rate is 
analogous to the difference in the direct doppler measured directly at the two 
sites with the use of coherent transponders on the spacecraft. In the doppler 
tracking technique, the phase change of the signal is continuously monitored 
over a full pass. The slope of the phase versus time plot yields the doppler at 
any particular time. The differential doppler reading is identical to the VLBI 
fringe-rate. The use of stable frequency standards in VLBI measurements 
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enables one to measure the differential doppler with improved resolution. 

However, it can not be ignored that the same information is obtainable from 
doppler tracking measurements made at the two terminals. The VLBI technique 
has certain advantages; (1) it is a passive technique, (2) an inexpensive wide- 
band noise source and a free running oscillator onboard the spacecraft are all 
that are needed, (3) ground station equipment is inexpensive and simple. With 
simplified data processing techniques to perform the cross-correlation, the 
VLBI measurements will become an operational tool for precision tracking and 
geodetic applications. 

One area of application we have in mind is to use the installations with 
large antennas to observe both satellites and natural radio sources. The quasar 
observations are used to calibrate the baseline and the station coordinates. The 
satellite VLBI observations of time-delay and fringe-rate are then combined with 
the calibrated baseline coordinates to obtain the spacecraft position very accurately. 
Then, small, portable terminals at remote sites are used to observe the much 
stronger signals from the spacecraft(s) [the satellite signals, in general, are 
six to seven orders of magnitude stronger than those from the celestial sources 
and permit high resolution time-delay interferometry] whose orbit(s) has been 
determined precisely relative to the intertial frame by VLBI observations from 
larger installations. The VLBI observations of satellite signals by the portable 
terminals are used to establish the geodetic ties directly. 

The philosophy outlined above was tried in the Rosman-Mojave experi- 
ment and later in the Rosman-Owens Valley- Agassiz experiment. The Rosman- 
Mojave baseline was determined with a precision of ±2.0 meters from VLBI 
observations of quasars. The calibrated baseline coordinates were then used 
in the ATS-3 satellite position determination. The Owens Valley- Agassiz base- 
line was determined with a precision of ±0.40 meters from quasar observations. 
(The quasar observations have been reported separately.) But no VLBI observa- 
tions have been conducted from remote sites using portable terminals. We will 
do so in our future experiments. The specific role of satellite VLBI measure- 
ments in the Earth and Ocean Dynamics studies will become clear in the very 
near future. 
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